Frequency multiplexed superconducting quantum interference device readout 
of large bolometer arrays for cosmic microwave background measurements 
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A technological milestone for experiments employing Transition Edge Sensor (TES) bolometers 
operating at sub-kelvin temperature is the deployment of detector arrays with lOOs-lOOOs of bolome- 
ters. One key technology for such arrays is readout multiplexing: the ability to read out many sensors 
simultaneously on the same set of wires. This paper describes a frequency-domain multiplexed read- 
out system which has been developed for and deployed on the APEX-SZ and South Pole Telescope 
millimeter wavelength receivers. In this system, the detector array is divided into modules of seven 
detectors, and each bolometer within the module is biased with a unique ~MHz sinusoidal carrier 
such that the individual bolometer signals are well separated in frequency space. The currents from 
all bolometers in a module are summed together and pre-amplified with Superconducting Quan- 
tum Interference Devices (SQUIDs) operating at 4K. Room-temperature electronics demodulate 
the carriers to recover the bolometer signals, which are digitized separately and stored to disk. This 
readout system contributes little noise relative to the detectors themselves, is remarkably insensi- 
tive to unwanted microphonic excitations, and provides a technology pathway to multiplexing larger 
numbers of sensors. 

PACS numbers: 07.20.Mc, 07.57.Kp, 95.55. Sh, 95.55. Rg 

Keywords: SQUID, multiplexed readout, bolometer, mm-wavelength astronomy 



I. INTRODUCTION 



A new generation of mm-wavelength receivers instru- 
mented with hundreds or thousands of transition edge 
sensor (TES) bolometric detectors, e.g., [IHSj, is allow- 
ing for unprecedented precision in measurements of the 
Cosmic Microwave Background (CMB). Large TES ar- 
rays are expected to bring forth similar advances in other 
observational bands such as the sub-mm (e.g., [6]), X- 
ray (e.g., [7]), and gamma-ray (e.g., [8HTn]) as well as 



phonon detection in dark matter direct detection exper- 
iments (e.g., jllj). 

TES detectors in the millimeter to far-IR range are 
approaching the photon statistics limit, meaning their 
noise performance is dominated by fluctuations in the 
arrival rate of photons at the detector, and not by the 
intrinsic noise performance of the detectors themselves. 
Consequently, higher mapping speed can best be real- 
ized by increasing the number of detectors, rather than 
by improving the detectors themselves. A key challenge 
in building focal planes with large sensor counts has been 
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multiplexing the signals between the cold sub-kelvin op- 
erating temperature of the detectors and room temper- 
ature electronics. Multiplexing minimizes the heat-load 
on the detector cold stage and reduces the complexity of 
the cold wiring. 

Two complementary multiplexing strategies are com- 
monly used: (1) time domain multiplexing (TDM) [121 - 
[T5] , where the detectors are read out sequentially one at 
a time with a revisit rate that is faster than the detec- 
tor time constant, and (2) frequency-domain multiplex- 
ing (fMUX) where each detector is biased at a unique 
location in frequency space and read out continuously. 
While the TDM system was developed earlier, both tech- 
nologies are now quite mature. An eight-channel TDM 
proof-of-concept system was first to make sky observa- 
tions [HI. A new TES readout scheme that multi- 
plexes signals from GHz frequency RF SQUIDs [171 El] 
is in the early stages of development, as is a new code- 
domain multiplexing readout [19| . 

In this paper, we describe the fMUX readout system 
[20l - f27] as it has been deployed for mm-wavelength ob- 
servations with the APEX-SZ 3 |28l [29] and South Pole 
Telescope (SPT) SZ [4, -'iO receivers. These experi- 
ments are designed for high angular resolution observa- 
tions of the CMB. For example, the unique combination 
of sensitivity and resolution of these instruments makes 
it possible to study the structure and evolution of galaxy 
clusters using the Sunyaev-Zel'dovich (SZ) effect [ST], a 
small distortion imprinted on the CMB when primor- 
dial photons inverse-Compton scatter off hot electrons in 
the intra-cluster medium. The resolved SZ image of 
the Bullet Cluster (IE 0657-56) obtained with APEX- 
SZ [32] was the first published scientific result with a 
large array of multiplexed TES detectors. SPT achieved 
the first SZ detections of previously unknown galaxy clus- 
ters [3S] with its array of multiplexed TES detectors. A 
wealth of new measurements have been presented by both 
SPT [SlHig and APEX-SZ 

Other implementations of frequency multiplexed TES 
readout systems are described in Refs. [50l I51j . 

This paper is outlined as follows. In ^Il| a brief 
overview of the system, its advantages, and the perfor- 
mance requirements is presented. A detailed descrip- 
tion of the sub-kelvin components of the system is pro- 
vided in §III| This description is accompanied by a sum- 
mary of the design process for specifying the basic mul- 



tiplexer parameters, such as channel spacing. { IV details 



the SQUID pre-amplifier and back-end electronics. The 
performance for the system operating in the SPT and 
APEX-SZ experiments is described in 3V| and conclu- 



sions are presented in S VI 



II. GENERAL OVERVIEW 

The fMUX readout system operates in concert with an 
array of TES detectors. Each bolometric detector has a 
metal absorber and TES, which are coupled by a weak 



thermal link to a 0.25 K thermal bath. A constant- 
amplitude sinusoidal voltage bias is applied to the TES, 
which is a superconducting film with a transition tem- 
perature of about 0.5 K. The combination of radiation 
power absorbed from the sky and electrical bias power 
raises the temperature of the TES mid-way into its su- 
perconducting transition. An increase in sky power 
incident on the absorber alters the TES resistance and 
(through negative electro-thermal feedback, see S III A I, 



produces a decrease in the electrical power being dissi- 
pated in the sensor. This, in turn, produces a change in 
the current through the sensor, since the applied voltage 
is held constant. TES thermometers are chosen because 
their steep transition allows for linearity, their monolithic 
wafer-scale fabrication yields large arrays, and because 
their low impedance makes them compatible with multi- 
plexing. The properties and parameters of the bolometer 
arrays dictate the design requirements for the readout 
system. 



A. The APEX-SZ and SPT Experiments 

The SPT focal plane [SH [5S] is divided into three 
separate bands centered at approximately 95, 150, and 
220 GHz, with bandwidths of about 40 GHz. The APEX- 
SZ detectors operate at 150 GHz with ~36 GHz band- 
width. Coherent amplification devices such as High Elec- 
tron Mobility Transistors (HEMTs) are ultimately lim- 
ited by quantum noise, and incoherent detection with 
bolometers can offer higher sensitivity at these frequen- 
cies. 

Both experiments use a closed-cycle refrigerator manu- 
factured by Chase Research Cryogenics (140 Manchester 
Road, Sheffield, UK SIO 5DL), with two ^He sorption 
stages and one ^He sorption stage, to cool the bolometers 
and sub-kelvin multiplexer components to their ~0.25 K 
bath temperature. The refrigerator is thermally coupled 
to the ^4K stage of a Cryomech (113 Falso Drive, Syra- 
cuse, USA 13211) model PT410 mechanical pulse tube 
cooler. With this combination, the system is cooled with- 
out the use of any expendable cryogens, greatly simpli- 
fying the operation of these instruments in the remote 
environments of the Chilean Atacama plateau and the 
geographic South Pole for APEX-SZ and SPT respec- 
tively. 

APEX-SZ and SPT have 0.4° and 1° fields of views 
(FoVs) respectively, and both experiments have full 
width half maximum (FWHM) beam sizes of about 1 ar- 
cminute. There are no "chopper mirrors" in the opti- 
cal systems to modulate the beams across the sky; in- 
stead the telescopes are scanned at typical speeds of 0.25- 
0.5°/second. Since we are interested in observing astro- 
nomical structures (such as primary CMB fluctuations) 
at scales occupying a large fraction of the FoV, these 
experimental parameters and scan strategy dictate that 
information will appear in the detector timestreams at 



frequencies as low as ~ 1 Hz. 
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FIG. 1: Layout of the frequency-domain multiplexed readout 
for a module of A'^ TES detectors. 



B. Readout System Concept 

The detector readout is segregated into identical fMUX 
modules, each biasing and reading out N TES detectors 
on a single pair of wires. The basic layout of one multi- 
plexer module of the fMUX system is shown in Figure [T] 
Figure[2]shows the system electronics (§ IV ) in greater de- 
tail. To set its specific bias frequency, each TES bolome- 
ter (-Rboio) is connected through a series-resonant circuit 
formed by an inductor (L) and capacitor (C) on the sub- 
kelvin stage. The inductance is the same for all chan- 
nels, defining the same i?boio/(27rL) electrical bandwidth 
for each detector. Thus, changing the capacitance 
alone sets the LCR^oio channel resonant frequency. A 
bank of N fixed amplitude sine-wave generators oper- 
ating at room temperature provides a 'comb' of detector 
voltage bias carriers, each operating at a frequency tuned 
to match the resonance of the individual LCRhoio chan- 
nels. Sky signals modulate the TES resistance, ampli- 
tude modulating the current of the carrier. This encodes 
the sky-signal in symmetrical sidebands above and below 
the carrier. Since the signal currents from the indi- 
vidual TES detectors are at different frequencies, they 
can all be summed together and transmitted to a sin- 
gle low-impedance input amplifier operating at 4K. We 
employ a series array SQUID operating in a flux-locked 
loop (ELL) for this transimpedance amplification stage. 
To reduce the SQUID dynamic range requirements, an 
inverted version of the carrier comb, referred to as the 
'nulling comb', is injected at the SQUID input to cancel 
the carriers (see ^ IV I . The comb of amplitude modulated 



carriers is transmitted from the SQUID output to a bank 
of room temperature demodulators. There is separate 
analog demodulator for each TES that mixes the detector 
signal down to baseband. A low-pass anti-aliasing filter is 
applied to the resultant timestream before it is digitized. 
The signal amplitude is proportional to the photon power 



deposited on the TES photon absorber. The APEX-SZ 
and SPT systems are configured with N = 7 detectors 
per readout module. 

The configuration shown in Figure [T] uses current sum- 
ming at the low impedance input of the transimpedance 
amplifier using a SQUID input stage. The fMUX can 
also be implemented by introducing a transformer in se- 
ries with each bolometer, and summing the voltages from 
the transformer secondaries |2I] , but since a TES biased 
by a constant voltage yields changes in signal current, the 
current summing configuration is more direct. The de- 
sign requirements, challenges, and implications are simi- 
lar for both configurations. For simplicity, we focus only 
on the current summing strategy. 

The fMUX readout architecture has the following key 
features: 

(1) Simplicity of cryogenic components: the sub-kelvin 
components consist of one capacitor and one supercon- 
ducting inductor per bolometer, followed by a single 
SQUID pre-amplifier per readout module at 4K. Since 
the SQUIDs are located far from the detectors, they can 
be magnetically shielded in a straightforward manner. 
The only custom sub-kelvin components are the induc- 
tors, which can be manufactured via a simple, low cost 
process with high yield. 

(2) Low wire count per detector and independent mod- 
ules. The fMUX requires just two wires per readout mod- 
ule for both bias and readout. Each module is completely 
independent so that wiring failures or other defects in one 
module do not affect others. 

(3) Signals are modulated at high (typically MHz) fre- 
quencies, well above microphonic pickup frequencies and 
amplifier low-frequency noise. The inductor-capacitor 
filters in series with the detectors suppress currents at 
frequencies outside the narrow LCRbolo bandwidth of 
the individual bolometer legs. This reduces crosstalk, 
eliminates the dominant source of microphonics usually 
present in bolometer systems, and results in a high rejec- 
tion of electromagnetic interference that comes from nu- 
merous sources such as scan-correlated changes in mag- 
netic field. 

(4) The readout system dissipates zero power on the de- 
tector cold stage as there are no resistive or active ele- 
ments. 

(5) Channels are read out continuously with no switching 
transients that limit the detector bandwidth. This makes 
the system easily applicable to higher bandwidth appli- 
cation such as photon counting X-ray microcalorimeters. 

(6) The voltage bias can be set to optimize the sensi- 
tivity of each individual detector. This allows the read- 
out system to compensate for detector fabrication non- 
uniformity. 

Practical challenges for implementing the fMUX sys- 
tem include: achieving the large dynamic range needed to 
accommodate the bias carriers; meeting the strict stabil- 
ity requirements on the bias generators needed to avoid 
introducing low frequency noise; managing non-idealities 
in cryogenic wiring that can introduce stray inductances 
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FIG. 2: Block diagram showing the basic topology of the fMUX readout system implementation, with an emphasis on the 
room temperature electronics. Refer to §IV| for a description. 



that spoil voltage bias or create cross-talk; avoiding de- 
tector oscillations which can affect other detectors in the 
same readout module through SQUID interactions; and, 
since the SQUID FLL includes a room temperature am- 
plifier, achieving loop stability in the present implemen- 
tation of this transimpedance amplifier is contingent on 
short wire lengths between 4K and room temperature. 



C. Design Requirements 

The TES detector properties, our telescope scan strat- 
egy, and our science goals define the following require- 
ments for the fMUX readout system: 
Noise: The noise contribution from the readout sys- 
tem should be small compared to photon and detector 
noise which is typically about 50 aW/\/Hz for ground- 
based mm- wavelength observations. The readout sys- 
tem receives noise contributions from its SQUID and 
room-temperature electronics amplifier stages, Johnson 
noise from the TES resistance, and Johnson noise leak- 
age through the LCRhoio fiher from neighboring detector 
channels. 

Detector Stability: The analog circuit that includes the 
TES must be stable across the full range of observation 
conditions. 

Cross-talk: Typical detector-to-detector cross-talk for 
our systems due to beam side-lobes and radiation leak- 
age in the bolometer integration cavity is ^ 1%. The 
readout system electrical cross-talk should be small in 



comparison. 

Dynamic Range: The readout system amplification needs 
to have sufficient dynamic range to accommodate the 
large bias carriers without introducing significant inter- 
modulation distortion. The dynamic range requirement 
is substantially reduced for the fMUX system by cancel- 
ing the carriers at the SQUID input with nulling signals 
that are inverted versions of the raw bias carriers, de- 



scribed in i IV 



Bandwidth and Post-Detection Sampling Rate: The tele- 
scope scan strategy defines a low frequency stability re- 
quirement of better than ~1 Hz which places stringent 
stability requirements on the sideband noise of the bias 
carriers. The fastest expected thermalization time con- 
stant of the detector absorber, r^eb —5 ms (often referred 
to as the "optical time constant"), defines 30 Hz as the 
upper end of the signal band. Note that much higher 
bandwidths are of interest for characterizing bolometers, 
since the response time of the TES thermistor can 
be much faster than the bolometer optical time con- 
stant. This motivates a sampling rate in the range of 
100-1000 Hz. Refer to TableUfor definitions of the time 
constants used in this paper and Ref. fSS] for a detailed 
thermal description of the detectors, including time re- 
sponse and diagrams. 

Programmahility: The system needs to operate au- 
tonomously under computer control, tuning and setting 
up each SQUID and detector. This tuning needs to be 
dynamic, so that changes in operating conditions (such as 
small changes in cryogenic temperatures or large changes 
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Twcb thermalization time constant of the detector radiation 
absorber 

TxEs effective TES response time for thermal signals 
tlcr electrical time constant of the notch filter formed by the 
bolometer resistance and multiplexer inductor/capacitor 

TABLE I: Definitions for the time constants utilized in this 
paper. 



in atmospheric loading) do not adversely affect the sys- 
tem yield or operation. 

Diagnostics: The system must provide sufficient diag- 
nostic information so as to characterize detector perfor- 
mance, recognize problems such as failed tunings, and 
provide feedback to the telescope operator. 




III. SUB-KELVIN MULTIPLEXER 

One of the advantages for the fMUX system is the sim- 
plicity of the sub-kelvin readout components, consisting 
only of capacitors, inductors, and wiring. However, as 
discussed in this section, successful implementation of an 
fMUX system requires careful management of stray re- 
actances to maintain good detector stability and optimal 
system noise. In addition to describing the specific im- 
plementations of the sub-kelvin multiplexer components 
for APEX-SZ and SPT in this section, after first describ- 
ing the TES bolometers, we outline the design process 
that leads to the choice of parameters. 



A. TES Bolometer Arrays 

The APEX-SZ and SPT TES bolometer arrays (Fig- 
ure [Sj) were fabricated at UC Berkeley [5iH5B] . The 
detector arrays consist of six triangular shaped wedges, 
which together form a hexagonal array of 330 and 966 
bolometers respectively for APEX-SZ and SPT. Short 
sections of circular wave-guide with conical feed-horns 
define the lower edge of the detector band and couple the 
detectors to free space. Metal mesh filters [57l|58] located 
above the feed-horns define the upper band edge. The 
bolometers are electrically connected to the cold mul- 
tiplexer components with aluminum wire bonds on the 
periphery of the bolometer array. 

The bolometers each consist of a 3 mm diameter gold 
absorber suspended above the 0.26 K cold stage bath 
through a thermal conductance in the range of G=100- 
200pW/K, where the bar indicates G is averaged across 
the temperature difference between the TES transition 
temperature and its thermal bath. The bolometer op- 
erates with total power Ptotai = ^(Tboio - ^bath) flow- 
ing to the bath. Optimizing G is a trade-off between 
(1) lower noise, which pushes G to lower values, and (2) 
high dynamic range for observations at a variety of op- 
tical loading conditions, which pushes G to higher val- 
ues. For APEX-SZ and SPT, the target G has been 



FIG. 3: Photograph showing one wedge of 161 SPT TES 
bolometers. The bolometer wedge is wire-bonded to a circuit 
board on which the multiplexer inductors and capacitors are 
mounted. An integrated fiexible microstrip (visible at the 
right of the photo) brings the leads to connectors mounted 
below the circuit board. The band-defining feed-horns and 
metal mesh filters have been removed for this photo. 



chosen so that the detector saturation power is slightly 
larger than twice the typical incident radiation power. 



total 



> 



2i^rad, which maintains high electro-thermal 
loop-gain (discussed below) under all observing condi- 
tions. 

The bolometer absorber has a spiderweb geometry, 
which provides a high optical efficiency while minimizing 
the cross section to cosmic rays. An Al/Ti bilayer TES is 
located at the center of each spiderweb and tuned to a su- 
perconducting transition temperature Tboio ^ 0.55 K and 
normal resistance i?normai 1.2 f2. The thermalization 
time of the spiderweb Tweb-^ typically 10-15 ms (16-11 Hz) 
for our detectors, sets the device's response time for op- 
tical signals. The TES has a response time r^gg to 
thermal signals that is sped up by electro-thermal feed- 
back and is typically more than an order of magnitude 
faster than Tweb- Since a fast T^f;g places additional 
constraints on the readout system (see S HI A 2 1 , the TES 



time constant is slowed down by attaching a gold feature 
with considerable heat capacity . Typical optical ef- 
ficiencies (for the conversion of power from astronomical 
radiation to electrical power) are ^80% for light incident 
on the detectors, and about 20-40% end-to-end from the 
telescope's primary dish. References [2] and [5^ provide 
discussions of the APEX-SZ and SPT optical cfficiences 
respectively. 



1. Detector Bias and Electro-Thermal Feedback 

The TES are biased with a fixed amplitude sinusoidal 
voltage at frequencies much higher (faster) than the TES 
can respond, so it sees the integrated bias power. The 
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steep Ai?boio/AT TES transition and voltage bias pro- 
vide the conditions for strong negative electro-thermal 
feedback (ETF) [SS]- When the radiation power P^ad 
incident on the bolometer increases, the bolometer tem- 
perature Tboio increases, driving up its resistance i?boio- 
This, in turn, brings down the electrical bias power 
-^eiec = ^1as/-^boio providcd by the sinusoid at frequency 
completing the negative ETF loop. 

ETF improves linearity, extends the dynamic range, 
and speeds up the effective TES response time for ther- 
mal signals. A measure of the ETF strength is the 
loop-gain, 



<-rJbolo 



(1) 



where a — (ilog(i?boio)/rflog(Tboio) is the logarithmic 
derivative of the bolometer resistance i?boio with tem- 
perature Tboioj G = rfPtotai/c^Tboio is the instanta- 
neous thermal conductance between the bolometer ab- 
sorber and bath, and Tboio is the TES transition temper- 
ature. When the loop-gain is high, the total power on 
the bolometer Ptotai — Prad + ^'^'cc "^i^^ ^e approximately 
constant and balanced by the rate of heat transmitted to 
the bath, Ptotai = G(rboio - Tbath) = Jrl^^°^° GdT . 

One simple strategy for providing voltage bias is to 
place the TES in a loop with a small bias resistor, 
^^bias << ^boio shown in Figure [l] across which the 
voltage bias is provided by an external current source. 
This loop consists of the bias resistor, one bolometer, the 
inductor and capacitor that define the loop resonant fre- 
quency, and the transimpedance amplifier input. Ideally, 
all other components in the loop would have negligible 
impedance in comparison with the bolometer. This 
includes the transimpedance amplifier that measures the 
TES current, wiring, and any other components such as 
the inductors and capacitors that define the resonant fre- 
quency. The wiring that connects the bolometers 
on the sub-kelvin stage to the low impedance ammeter 
(in this implementation, a SQUID array located on a 
4K stage) can be substantial in length and span large 
temperature gradients. Wiring resistance can be kept 
low while maintaining low thermal-conductivity by con- 
structing the wires from a superconductor, such as NbTi. 
Since this wiring lies outside the tuned LCPboio leg of the 
loop, its inductance L can not be tuned out by the fMUX 
capacitor for more than one channel. This means it is 
necessary to keep the wiring inductance low enough that 
the reactance is small compared to Pboio at the resonant 
frequency. Meeting this requirement is difficult with tra- 
ditional twisted pair wiring (50 nH amounts to O.SjTi at 
IMHz, where j denotes the imaginary component), so 
that low inductance microstrips are necessary for all or 



part of this wiring (see SIVD). 



8. Bandwidth and Detector Stability 

To avoid under-damped or exponentially growing os- 
cillations in the TES response, the electrical power (half- 
width at half maximum, HWHM) bandwidth of the 
iCPboio filter 1/tlcr must exceed the TES thermal 
bandwidth l/rf^s by a factor I/tlcr > 5.8/T^%a, fSU] 
(see Ref. [61] for a thorough description). For the 
APEX-SZ and SPT detectors with Pboio ~ 1 and 
thermalisation time constant r^cb > 5 ms (32 Hz roll- 
off frequency), we employ 15.8 /iH inductors to maxi- 
mize the detector readout bandwidth while still meeting 
our constraint on cross-channel current leakage. This re- 
sults in the requirement r-f-gg > 0.2 ms (0.8 kHz roll-off 
frequency). This allows for an order of magnitude range 
(0.2 ms < r^Eg < 5 ms) for the TES time constant. The 
TES time constants are increased to lie within this range 
by adding a heat capacity in the form of an Au island 
that has a strong thermal coupling to the TES. How- 
ever, since the t^%s is sped up by the loop-gain of the 
ETF, it is important to have a substantial stability range. 
This stability requirement defines the minimum allowed 
LCRhoio bandwidth, and hence the channel spacing (de- 
scribed below) for the fMUX system. 



B. Multiplexer Channel Spacing 

The spacing between adjacent multiplexer channels is 
chosen such that (1) Johnson noise leakage from adjacent 
channels is small, (2) bias current leakage from adjacent 
channels does not spoil the detector voltage bias, and 
(3) cross-talk between neighboring detectors on a comb 
is small in comparison to the other sources of cross-talk 
in the experiment. These conditions must hold for the 
LCRhoio bandwidth that is defined by the ETF stability 
requirements discussed in §111 A 2| For APEX-SZ and 
SPT a channel spacing of 75 kHz is used. The first two 
criteria are discussed below, while cross-talk is discussed 
in j |HTC| 

Each bolometer in the multiplexer module is a source of 
broadband Johnson noise. The SQUID input receives the 
sum of all Johnson noise contributions from every chan- 
nel, so it is necessary to space the channels far enough 
apart in frequency so that the Johnson noise from neigh- 
boring channels is attenuated to a negligible level. This 
requirement is not the dominant constraint for channel 
spacing. At a frequency of 75 kHz away from resonance, 
for L = 15.8 /xH and Pboio = the Johnson noise 
current of a detector is attenuated by a factor 15, pro- 
viding an increase of just 0.2% to the total Johnson cur- 
rent noise of its neighbor. The increase in total detector 
noise is substantially smaller than this when other noise 
sources are also considered. The Johnson noise leakage 
from spacings as close as 4 times the LCPboio bandwidth 
(20 kHz, in this case) or less could be tolerated, con- 
tributing a 3% increase in the Johnson noise current of a 
neighbor. 
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A more serious constraint on channel spacing is bias 
current leakage. In addition to its own on-resonance 
bias, bolometer channel i also sees a fraction of the cur- 
rent from the i ± 1 nearest carriers that neighbor in fre- 
quency space. The magnitude ratio of the off-resonance 
current IqI^I from the neighboring carrier bias to the on- 
resonance current /^jj,; flowing through bolometer chan- 
nel i is 



'Chi 

TUli 

-'Chi 



-Rbolo 



(2) 



where i?boio is the bolometer resistance, L is the series in- 
ductor, Cchi is the series capacitor, and uj.i±i is the angu- 
lar frequency of the adjacent bias carrier. For the channel 
spacing and parameters considered here, this may be ap- 
proximated IqI^I l^ch i ~ ^Sl ' 'wh^re Aw is the channel 
spacing. 

Bias carrier leakage can be problematic, mostly be- 
cause of cross-talk (discussed below), but also because 
the leakage is a quasi-current bias that provides a pos- 
itive ETF contribution to the total bias power seen by 
the bolometer. While the negative ETF afforded by the 
on-resonance voltage bias will dominate, the contribu- 
tion from carrier leakage can cause instability and will 
also enhance the bolometer responsivity slightly. 

All of these effects are rather small in practice. For a 
spacing of 50 kHz (10 times the LCi?boio bandwidth), the 
electrical bias power from one of the nearest neighbors 
is roughly 1% of the on-resonance power assuming that 
all bolometers are biased at the same voltage amplitude. 
For channels biased near the center of the bias-frequency 
comb, the total off-resonance bias power is approximately 
3%. This leads to a 3% increase in both the responsitivity 
and the effective time constant, and a 3% decrease in the 
maximum loopgain. At our spacing of 75 kHz (15 times 
the LCBhoio bandwidth), these effects are reduced to 1%. 



C. Cross-talk 

The multiplexer contributes channel-to-channel cross- 
talk due to (1) inductor cross-coupling, (2) bias carrier 
leakage, and (3) signals from a detector causing heat- 
ing in neighbor detectors due to the voltage drop across 
non-zero SQUID or wiring impedance. Contributions 
from these mechanisms, discussed below, are smaller than 
the ^1% optical cross-talk that typically exists between 
neighboring pixels for APEX-SZ and SPT. 

(1) Inductor cross-coupling: The mutual inductance 



ki,jLiLj between inductors Li and Lj of two 



channels i and j in a multiplexer module creates cross- 
talk. The carrier current li flowing through the i*'* 
bolometer channel induces a voltage \Vj\ = uiiMijIi in 
the inductor j that physically neighbors channel i. As 
changes in response to sky signals, the voltage Vj 
is modulated as well. This modulation produces small 
changes in the current flowing in channel j since Rl,^^^ 



forms a small fraction of the total LCR^^^^^ impedance 
at frequencies far from its resonance. The measured cou- 
pling coefficient |25j between neighboring inductors on 
the same chip is A: = 0.010 ± 0.002 for the fMUX system. 
This form of cross-talk can be made negligible by ensur- 
ing that channels with physically neighboring inductors 
are not neighbors in frequency space (i.e., that physically 
neighboring inductors in this implementation have carrier 
bias frequencies separated by more than 150 kHz). 

(2) Bias carrier leakage: As i?boio changes in response 
to sky signals incident on bolometer «, both its on- 
resonance carrier at frequency uji and the off-resonance 
leakage bias carriers at neighboring frequencies uii±i are 
modulated, creating cross-talk. This mechanism encodes 
a cross-talk sky-signal on the leakage current waveform. 
The signal appears as an amplitude modulation in the 
on-resonance carrier, but leads to modulation in both 
the amplitude and the phase of the off-resonance carrier. 
Thus a full quadrature demodulator would see cross talk 
signal in both the I- and Q-components of the demodula- 
tor output. However, the demodulator measures only the 
I-component modulations. For an off-resonance LC-Rboio 
channel i ± 1 at carrier frequency oJi , the amplitude fluc- 
tuation is: 



bias 



-Rbolo + i^^iL + l/(jWjCchi±l) 



-Rbolo + j2AujL 

jRholo 



bias 

j2ALuL 



1 



2AwL 



(3) 
(4) 
(5) 



where a first-order Taylor expansion in the parameter 
has been performed for the last equality. This re- 



2Aw L 



suits in a current modulation with a change in bolometer 
resistance of 



A/; 



Chi±l 



Ai?boio (2AwL)2' 
which should be compared to the current modulation 



(6) 



AT. 



Chi 



bias 



Ai?bol( 



^bolo 



(7) 



that would occur on-resonance. The magnitude ratio of 
Equation [6] to Equation [7] 



n2 

^bolo 



(2Aa;L)2 



(8) 



is a good approximation for this cross-talk. This is 
equivalent to the square of the current leakage ratio 
presented in Equation [2] For our circuit parameters 
(i?boio = 0.7517, Aw = 27r • 75 kHz, L ^ 15.8 //H), this 
effect results in ^0.25% cross-talk. Note that this is 
larger than the cross-talk estimate that was derived 
in Ref. [25,, which did not correctly take into account 
the phase of this cross-talk signal. 
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(3) Non-zero SQUID/wiring impedance: A change in 
bolometer impedance of channel i in response to sky- 
signals, causes the voltage at carrier frequency uJi across 
the stray SQUID and wiring reactance to be modu- 
lated, which in turn causes a modulation in the leakage 
current power deposited in neighboring detectors i ± 1. 
This mechanism transmits cross-talk power from the on- 
resonance detector to the off-resonance neighbors, which 
is the opposite transfer direction compared to mecha- 
nism (2) . The voltage V,noduic across the sub-kelvin mul- 
tiplexer module (measured from the right side of pi to 
the right side of p2 in Figure |4| is not exactly equal to 
the voltage Vbias across the bias resistor. This discrep- 
ancy arises from the non-zero complex input impedance 
of the SQUID system (p5 in Figure [4]) and the stray in- 
ductance of the wiring (pi and p2 in Figure |4]) that con- 
nects the module to the SQUID input (together labelled 

14tray) SUCh that Kiodulc = ^bias " 14tray For the pur- 

poses of this discussion, the p3 and p4 stray components 
in Figure [4] can be ignored. Consider the current /^j^^ 
at carrier frequency that flows predominantly through 
the i*'' bolometer. The change in voltage across the mod- 
ule due to a change d/^h i 



dVr 



module 



c'^bias — c'Kitray 



^C^Kitray — ^'^-^Ch i-^stray 



-dP, 



Chi 



stray 



(9) 
(10) 

(11) 



at fixed Vbias- We have assumed that the SQUID sys- 
tem and wiring impedances Zgtray are dominated by an 
inductive term Lgtray This voltage induces a leakage 
current dl'^-^^^-^ through the neighboring i± 1 bolometer 
channel 



dl[ 



Chi±l 



'^I'niodulc dl^^j^ ■ jUJiL^fcay 



7LCR 



j2AwL 



(12) 



where 2'^hi±i — j'^^'^L is the impedance of the neigh- 
boring i ± 1 LC Rhoio leg of the cold multiplexer module 
at frequency uJi. This current fluctuation deposits Joule 
heating power 
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FIG. 4: Stray reactances in the sub-kelvin multiplexer circuit 
are problematic since the circuit operates at MHz frequency 
with characteristic impedances ~ Ifi. The stray inductance 
of the wiring connecting the bias resistor (-Rbias) to the mul- 
tiplexer circuit is shown at pi and p2, the stray inductance 
in series with -Rbias is shown at p3, the stray capacitance in 
parallel with a channel's inductor is shown at p4, and the non- 
zero input impedance of the transimpedance pre-amplifier is 
shown at p5. 



While the three forms of readout system cross-talk dis- 
cussed above would not be present in a non-multiplexed 
system, we note that a dominant form of electronic cross- 
talk for a non-multiplexed system can be eliminated for 
the fMUX system. Capacitive coupling of wires within 
multi-wire cables that carry signals from different detec- 
tor channels is typically ~1% for bundles of twisted pairs. 
For a non-multiplexed system, all detector signals nor- 
mally occupy the same frequency bandwidth as they are 
transmitted through these cables, resulting in 1% cross 
talk with neighbors. For the fMUX system, each wire 
carries a module of detectors, each at a unique bias fre- 
quency. Since the ^100 Hz signal bandwidth is a small 
fraction of the total bandwidth, the bias frequencies of 
modules that are bundled together can be chosen to be 
unique. This places the cable cross-talk signals in a por- 
tion of the bandwidth that is discarded by the signal de- 
modulation. This means that, overall, this multiplexed 
readout system has lower cross-talk than a typical non- 
multiplexed system. 



d'Pchi±l — 2-Rbolo -^Clii±l ' ^^Cli j±l (13) 

in the neighboring bolometer channel, which should be 
compared to the signal in Channel i, 

dPchi — dichi Kiodulc- (14) 
The ratio of these power fluctuations is the cross talk, 

'^Pchi±l -^Ch'i±l "^i -^stray 



dP'^' 



(15) 



For the fMUX system with L/Lgtray ^ 150, -Rboio = 
0.7517, Aw = 27r - 75 kHz, and taking = 27r - 750 kHz, 
this cross-talk amounts to 0.3%. 



D. Channel Count 

Considering the results of the previous sections, the 
number of channels that can be multiplexed together 
without adversely affecting the noise or cross-talk for a 
system with a given set of parameters is derived as fol- 
lows. 

The LCRhoio HWHM power bandwidth is chosen to 
be at least 5.8 times wider than the inverse of th e fastes t 
expected TES time constant, l/(27rr^|g) (see § |lIIA2| ). 
This has very conservatively been chosen to be 5 kHz 
in this system, allowing an order of magnitude range of 



cff 
TES 



0.2 ms < T, 

ordered appropriately as 5.8tlcr < t-tes ''wob 



< 5 ms so that the time constants are 

cff 
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The spacing between LCRhoio resonances is chosen so 
that the total bias carrier leakage current is less than 20% 
of the on- resonance carrier current, and the cross-talk is 
less than 0.5%. For the parameters of this system, this 
results in a spacing of 15 times the LCRhoio bandwidth, 
75 kHz. This spacing also keeps Johnson noise leakage 
at a negligible level. 

Having specified the channel spacing, the number 
of channels is defined by the bandwidth of the low 
impedance amplifier. The SQUID shunt-feedback cir- 



cuit (discussed in [ IV B I for this system has a bandwidth 
of 1 MHz. The carriers are positioned between about 
400 kHz and 900 kHz, ensuring odd harmonics from low- 
frequency carriers do not fall within the bandwidth of 
higher frequency channels. Odd harmonics are the dom- 
inant form of distortion from SQUID non-linearity. Bias 
carriers below 400 kHz are avoided because of the diffi- 
culty of finding high-Q capacitors with sufficiently high 
capacitance that behave well at cryogenic temperatures. 
These considerations result in 7 channels per multiplexer 
module so that the number of wires connecting the sub- 
kelvin multiplexer to the 4 K SQUID is the total number 
of detectors multiplied by 2/7. 

The per-module channel count for APEX-SZ and SPT 
is conservative. A small increase in channel count could 
be achieved by slowing the TES time constants t^%s so 
that they are better matched to the absorber thermaliza- 
tion time r^ch, allowing narrower channel spacing. Sub- 
stantial increases could be achieved by implementing a 
SQUID fiux-locked loop (FLL) with substantially larger 
bandwidth such as Ref. |62j . Finally, much narrower 
channel spacing could be realized by using a sharper 
band-pass filter than the single pole LCRhoio employed 
here, although this would substantially complicate the 
cold circuit. 



E. Sub-kelvin Multiplexer Implementation 

Each sub-kelvin multiplexer module consists of a 
custom-fabricated chip of seven inductors, seven commer- 
cially available ceramic capacitors with negative-positive- 
zero (NPO) dielectric material, and a custom printed cir- 
cuit board (shared by all multiplexer modules for a par- 
ticular detector wafer) on which these components are 
mounted. The multiplexer modules are interfaced with 
aluminum wire bonds to the detector wafer as shown in 
Figure [3j The circuit board is thermally sunk to the same 
250 mK temperature stage as the bolometers. 

The inductors were custom fabricated at the Northrop- 
Grumman superconducting micro-fabrication facility 
with a lithographic process. Each chip consists of eight 
inductors, of which we use seven. One inductor is about 
1.5 mm X 1.5 mm, consisting of a 140 turn spiral (2/im 
width niobium trace with pitch of 4 fiTo) separated by a 
50 nm Si02 insulating layer from a niobium, flux-focusing 
washer that is left floating. The 15.8 /iH inductance is 
the same for all channels in the module. The Northrop- 
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FIG. 5: The resonant frequencies for each LCRhoio leg of 
the multiplexer circuit are measured by frequency-sweeping 
a small amplitude bias carrier across the range of carrier 
frequencies with the bolometer stage temperature held just 
above the detector superconducting transition. The resultant 
network analysis data, shown in the amplitude vs. frequency 
plot above, is fitted with a model to determine the optimum 
carrier bias frequency for each bolometer. One such fit is 
shown in the inset. 



Grumman superconducting microlab is no longer avail- 
able, and inductors for new implementations of the fMUX 
system are now being fabricated at the National Institute 
for Standards and Technology (NIST). 

The frequency of each channel resonance is set by the 
series capacitor, for which we use commercially available 
surface mount capacitors. Our target spacing for the 
7-channel multiplexer is 75 kHz. Several capacitors, pre- 
selected by measuring their capacitances individually, are 
stacked in parallel to achieve the target capacitance for 
each channel. The capacitances change by several per- 
cent when cooled from room temperature to 0.25 K, but 
in general all devices shift proportionately and the chan- 
nel spacing is sufficiently preserved. 

The resonant frequencies for each multiplexer module 
are measured by frequency-sweeping a small amplitude 
bias carrier across the range of carrier frequencies with 
the bolometer stage temperature held just above the de- 
tector superconducting transition. The amplitude of the 
current through the SQUID is measured. A representa- 
tive 'network analysis' amplitude vs. frequency plot for a 
multiplexer module is shown in Figure [5] Each peak rep- 
resents one bolometer channel, with its location, width, 
and height defined by the LCRhoio parameters. The bias 
frequency for each bolometer is determined by fitting a 
model to each peak that includes parameters to account 
for the current flowing through off-resonant LCRhoio legs 
in the multiplexer module. 
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F. Effects of Stray Reactances 

Although the sub-kelvin muhiplexer components are 
relatively simple, non-idealities in the form of stray in- 
ductances and capacitances can be problematic, since the 
circuit operates at MHz frequencies with characteristic 
impedances of about an ohm. The effect of these strays 
are discussed below. 

The inductance of the cryogenic wiring that connects 
the bias resistor (i?bias) at 4 K to the sub-kelvin LCRhoXo 
multiplexer circuit is the most important stray, shown at 
pi and p2 in Figure |4j This inductance is outside each 
LCi?boio resonant circuit, and so cannot be tuned out 
at all of the carrier resonant frequencies. The voltage 
drop across this reactance reduces the voltage across the 
bolometer and partially spoils the constant voltage bias 
that is necessary for electro-thermal feedback. The ef- 
fect is worse at higher frequency. To control this stray 
component, the wiring is constructed of two segments: a 
low inductance (^0.3nH/cm) broadside coupled stripline 
is used for the majority of the wiring distance, followed 
by a short segment of NbTi twisted pair (~7nH/cm) 
that is necessary to create a low thermal conductiv- 
ity gap between the two temperature stages. The to- 
tal inductance for this wiring in the system implementa- 
tion is about 105 nH, amounting to (0.2-0.6) jTi at 300- 
900 kHz, increasing the impedance of the on-resonance 
0.75 bolometer circuit (-Rboio + j^Lgtray) by 3-27%. 

There is also a much smaller stray inductance in series 
with i?bias (p3 in Figure [i]) that has the opposite effect 
compared to the above, increasing the bias voltage across 
the bolometer. The circuit board housing the SQUIDs 
and 30 mf2 bias resistors has been optimized to minimize 
this stray, which is estimated at ~4 nH. The total Rbias + 
jojistray rcactauce increases from 31 voTl at 300 kHz to 
SSmfi at 900 kHz. While this stray results in a large 
change (as much as 26% at 900 kHz) in the effective bias 
voltage across the bolometer, this phase and amplitude 
modified voltage still provides an effective bias for the 
purposes of ETF. To see this, we note that it would be 
possible to bias the bolometer using just an inductance 
in place of the bias resistor, eliminating Johnson noise 
from this component. The drawback is that the strong 
dependence on frequency of this inductance-bias voltage 
strategy would complicate the setup and tuning of the 
system. An inductive or capacitive voltage divider, as 
described in Rcf. 50J, would address both issues. 

Both stray inductances described above affect the 
phase of the bias carriers in addition to their amplitude. 
Fortunately, the phase shifts of the two effects partially 
cancel. However, because of the lack of phase tuning ca- 
pability for the demodulator in the backend electronics 



in this implementation of the readout (S IV C I, this phase 



shift contributes to a slight degradation of noise at high 
frequency. 

A stray capacitance is in parallel with the inductor (p4 
in Figure |4]) that forms the LCi?boio series resonance for 
each multiplexer channel. This is due to a flux-focusing 



washer below the spiral trace of each inductor. With 
the present architecture, the resonance is in the 16 MHz 
range. While this feature is well above the multiplexer 
bandwidth, the phase shift it produces can cause stabil- 
ity problems for the SQUID pre-amplifier feedback, dis- 
cussed in |IVB[ 

The transimpedance pre-amplifier (a SQUID array in 
this implementation) has non-zero input reactance (p5 in 
Figure [4|. This means that currents injected along the 
nulling wire are split between two return paths: (1) the 
desired path through the pre-amplifier, and (2) through 
the bolometers. We refer to path (2) as "nuUer leakage" . 
Like the carrier voltage across the bias resistor, this leak- 
age also appears as a voltage across the bolometers, with 
the pre-amplifier reactance playing the role of a bias re- 
sistor. This effect is more pronounced at higher frequen- 
cies. When the nuller current phase and amplitude is ad- 
justed so that the net (carrier and nuller) current through 
the pre-ampifier is zero, the nuller current will be an ex- 
act inverted copy of the current through the bolometer. 
The nuller current is adjusted in this manner when the 
bolometers are tuned, forming a 'virtual ground' at the 
pre-amplifier input so that the net voltage across the in- 
put is zero after each tuning. Changes in observation 
loading conditions between bolometer tunings will result 
in a non-zero voltage across the pre-amplifier and alter 
the voltage bias across the bolometer. 

Achieving the design requirements outlined in § |II C| 
has required careful handling of the stray components 
described above. Future readout systems that employ 
higher frequency carriers to multiplex a larger number of 
bolometer channels will need more advanced techniques 
to correct for these stray components. 



IV. SQUID PRE-AMPLIFIER AND ROOM 
TEMPERATURE ELECTRONICS 



Having described the design considerations for the sub- 
kelvin components in the previous section, we now de- 
scribe the SQUID control and backend electronics in this 
section, beginning with the generation of the carriers 
that bias the detectors, then describing the SQUID pre- 
amplifier and its room temperature control electronics, 
and ending with the demodulation electronics. 

The topology of the system is shown in Figure [2] 
for one 7-channel fMUX module. APEX-SZ uses 40 
fMUX modules and SPT uses 120. Seven bias carri- 
ers are synthesized separately using Direct Digital Syn- 
thesizers (DDS) located on room temperature oscil- 
lator/demodulator boards. These analog carriers are 
added together to form a bias 'comb' and transmitted via 
a twisted pair cable and the room temperature SQUID 
controller board to the cryostat, where a 30 mO bias re- 
sistor at 4 K produces a voltage bias across the bolome- 
ters operating at 250 mK. LC resonant filters in se- 
ries with the bolometers select a single carrier tone for 
each bolometer. Changes in optical power on the de- 
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FIG. 6: (Color online) Photo showing the oscilla- 

tor/demodulator circuit board that provides the biases and 
demodulates the outputs for 14 bolometer channels. 



tectors induce amplitude modulation of the carrier cur- 
rents which are summed together and pass through the 
input of a series-array SQUID device. The SQUID is 
coupled to a room temperature operational amplifier to 
form a shunt-feedback flux- locked- loop circuit. A sep- 
arate 'nulling comb', which is simply an inverted ver- 
sion of the carrier comb, is transmitted from the oscilla- 
tor/demodulator boards to the SQUID input. This serves 
to cancel the large carriers, reducing the dynamic range 
requirement for the SQUID. The nulled and amplified 
comb is transmitted via a cable to a bank of seven demod- 
ulators located on the oscillator/demodulator boards. 
Each amplitude-modulated carrier is demodulated sep- 
arately, low-pass filtered, and digitized. The digitized 
sky-signals are transmitted to a data acquisition com- 
puter and recorded on disk. Each of these components is 
described in detail below. 



A. Carrier Generation 

The bolometer bias carriers are synthesized on custom 
circuit boards (Figure |6| that each handle 14 bolometer 
channels. In addition to synthesizing the carriers, these 
boards also perform the demodulation and digitization of 



the bolometer outputs described in { IV C The layout of 
these boards is optimized to minimize cross-talk between 
channels and pickup from other electronics systems. A 
block diagram showing the circuit used to generate the 
carrier comb and nuUer comb is shown in Figure [Tj 

The bias for each bolometer is generated with a Di- 
rect Digital Synthesizer (DDS, Analog Devices AD9854) 
that produces a sine wave by storing the waveform in- 
stantaneous phase in an accumulator and incrementing 
this accumulator every clock cycle. It uses the phase 
to address a corresponding amplitude in a look-up ta- 
ble, and sends that amplitude to a digital-to-analog con- 
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FIG. 7: Block diagram showing the carrier and nuller gener- 
ation circuits. The sine and cosine outputs at the DDS have 
amplitude adjustment registers that are not shown. The DDS 
amplitudes are used to control the phase of the nuller signals 
by adjusting the admixtures of sine and cosine waves. 



vertcr (DAC) sampling at 40 MHz. Because the phase 
information is held in only one register, it is not costly 
to achieve many bits of frequency resolution (these de- 
vices use a 48-bit phase accumulator and have fiKz reso- 
lution). Having accurate amplitude information is more 
costly since each memory location within the look-up ta- 
ble needs this accuracy. Our devices provide 12-bit am- 
plitude information. The output waveform can be multi- 
plied pre-digitization by an amplitude adjustment factor 
that ranges from to 1. In addition to the DAC provid- 
ing the sine wave output, each DDS has a second DAC 
output that has a fixed 90 degree phase off'set (i.e., a 
cosine wave) and separate amplitude adjustment factor. 
The DDS chips are the single biggest power consumers 
in the readout system, each one consuming ^0.4 W. 

The DDS output signals are differential and each pass 
through low-pass 5.5 MHz analog filters to reduce clock 
bleed-through and limit the noise bandwidth of the bias 
carriers. The sine-wave output is used for the bolome- 
ter bias signals. An 8-bit digital potentiometer is used 
to control this amplitude. A coarse voltage attenuator 
can be enabled to provide further amplitude adjustment, 
allowing for good noise performance over a broad range 
of bolometer parameters. The analog differential out- 
put currents from eight oscillator circuits are summed 
together to produce a bias comb that is transmitted to 
the cryostat via a twisted pair cable and the SQUID con- 
troller boards. Seven of these eight outputs are used to 
bias the bolometers, the eighth is used as a test signal 
that can be injected one at a time to each of the bolome- 
ters in the comb. 

The full scale carrier to white noise ratio is « 6 x 
10^\/Hz, dominated by the digitization noise of the DDS. 
More important than the white noise level however, is the 
carrier sideband noise, since this noise eventually domi- 
nates at low frequencies and limits our ability to observe 
large angular scale signals on the sky. Carrier sideband 
noise can originate from clock jitter, DDS voltage refer- 
ence jitter, and/or current noise in the transistors used 
in the DDS DAC. This system is not sensitive to clock- 
jitter, because the demodulator is driven by the same 
oscillator as the bias, approximately canceling the effect. 
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The dominant source of sideband noise is the low fre- 
quency noise of the DAC output transistors. This noise 
is proportional to the amplitude Ac of the carrier and 
measured to be ^ x 10~^/\/Hz at 1 Hz away from the 
carrier, with a 1/f power spectrum. 

In addition to the bias carrier generation, the same 
DDS circuits are used to provide the nulling signal that 
zeros the carrier amplitude just before the SQUID input. 
The nuUer consists of an inverted version of the bias sine- 
wave, with a small admixture of the cosine signal from 
the second DDS output. The sum of these two signals is 
also a sine wave with an phase that can be programmed 
(by adjusting the cosine amplitude) from about 158 to 
202 degrees from the bias carrier. This phase adjust- 
ment compensates for analog phase shifts between the 
carrier and nuUer that arise because the two signals take 
different paths through the cold multiplexer (Figure |2]). 
The filtering, amplitude adjustment, coarse attenuation, 
and summing circuits are the same for the nulling signals 
as for the carriers. 

This scheme for synthesizing the nulling signal was mo- 
tivated so that, when the cosine admixture is small, the 
bias low-frequency carrier sideband noise would be in- 
verted for the nuller signal and subtract to zero at the 
input of the SQUID. This sideband-noise nulling was 
demonstrated to work very well when adding the carrier 
to the nuller through a resistor. However, the scheme 
does not work when used to bias bolometers, since the 
bolometer has an effective negative resistance within its 
active bandwidth. This means the sideband noise of the 
carrier is inverted by the bolometer, while its MHz car- 
rier is not inverted. The net result is that the carrier 
sideband noise adds coherently to the nuller sideband 
noise. Though this results in excess low- frequency read- 
out noise for this system, other sources of low-frequency 
noise dominate our detector noise (e.g., residual atmo- 
sphere, temperature drifts, etc.), and the low frequency 
noise from the readout does not dominate the overall low 
frequency noise of the experiment. The digital fMUX 
system described in Ref. [5T]) uses a different scheme for 
nulling that avoids this excess readout noise. 



B. SQUID Electronics 

Transimpedance amplification of the bolometer cur- 
rents is achieved with a series-array SQUID operating 
in a shunt-feedback circuit with a room temperature op- 
amp as shown in Figure |8] This circuit has the following 
properties: (1) Its input impedance is sufficiently low so 
as not to spoil the voltage bias across the bolometers. (2) 
It provides sufficient transimpedance so that the output 
signals can be interfaced to standard room-temperature 
amplifiers. (3) Its noise, referred to the input, is small 
compared to the bolometer noise. (4) It is sufficiently 
linear that intermodulation distortion is not an issue for 
the system operation. 

For the fMUX system, the SQUIDs are housed in 



4 K 



Room 

Temperature 




FIG. 8: The SQUID shunt-feedback circuit is shown. The 
capacitor and resistor that together form the lead-lag filter 
(refer to text) connect the non-inverting input of the opera- 
tional amplifier to ground. 



groups of eight on custom SQUID-mounting boards that 
are heat-sunk to the 4 K cryostat mainplate. The rest of 
the shunt-feedback circuit resides on custom room tem- 
perature SQUID-controUer circuit boards, each of which 
manages eight SQUIDs. In addition to completing the 
shunt-feedback circuit, these electronics provide the bias 
currents for the SQUIDs and condition all analog signals 
entering the cryostat by attenuating them and filtering 
them for RF. 



1. Series Array SQUID 

Each bolometer module is connected directly to the 
input coil of a series- array SQUID manufactured at 
NIST Each array consists of 100 individual SQUID 
elements with their input coils, feedback coils, and out- 
puts connected in series. The input signals add construc- 
tively while the output noise voltage of the individual 
elements adds incoherently, resulting in a signal-to-noise 
ratio that is enhanced by y/n, where n = 100 is the num- 
ber of elements. The array noise is 3.5 pA/vT5z referred 
to the input coil and the bandwidth is 120 MHz. 

Each SQUID element has two inductor coils that were 
intended for use separately as input coil and feedback 
coil. For shunt-feedback, the input coil is used for both 
functions. This coil array has an input inductance Lj^ ~ 
150 nH and mutual inductance for each 8-turn element 
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The response function for an open-loop SQUID may 
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be approximated as sinusoidal, with output voltage 

V^l^j. ~ VpK sin(27r//Ar//$J, (16) 

where Ijn is the current through the input coil, Vpx is 
the peak output voltage of the device, and /$o is the input 
current required to provide a quantum of flux $o = h/2e 
through the SQUID. For our devices /$q « 25 fj,A and 
VpK ~ 1.5 — 3.5 mV depending on fabrication repro- 
ducibility, operating temperature, and tuning. One can 
choose the operating point along the sinusoidal response 
curve using a static flux bias. For the shunt-feedback 
circuit, the devices are operated on the mid-point of the 



falling edge (see ^ IV B 5 1 . The derivative of Equation \16\ 
at this point is the transimpedance, Z^^^^^ ~ —500V/A, 
in good agreement with the measured performance of the 
devices. 



2. Shunt-feedback Circuit 

Shunt-feedback (Figure [s] ) is essential for linearizing 
the SQUID response and reducing the input reactance of 
the FLL circuit so that the bolometer impedance dom- 
inates the bias loop. Had the more common feedback 
configuration with separate feedback and input coils been 
used, it would increase the FLL circuit input impedance 
substantially, spoiling the bolometer voltage bias and rel- 
egating it to an unstable regime of electro-thermal feed- 
back. 

Shunt-feedback is rarely used for practical SQUID ap- 
plications, because it couples the the strong input coil di- 
rectly to the room temperature electronics, where pickup 
from digital components or RF signals are more difficult 
to shield against. Implementation of a shunt-feedback 
system requires careful filtering and shielding to prevent 
outside unwanted signals from reaching the SQUID input 
coil through the feedback line. Since this implementation 
is uncommon, the basic equations that characterize the 
shunt-feedback circuit operation are presented below. 

An open-loop SQUID has input reactance Zf^ ~ 
jujL^^ (the reflected impedance of the SQUID is small 
at these frequencies and has been omitted), where uj is 
the angular frequency. The forward voltage gain of the 
open- loop SQUID/amplifier cascade is 
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where Vq^^ and Vjff are the voltages at the ampli- 
fier output and across the SQUID input coil respec- 
tively, Ij^ is the current through the SQUID input 
coil, and G^^^ ~ 250 is the non-inverting amplifier 
forward gain. The SQUID transimpedance is written 
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trans — ~\^trans\ make explicit its operation in in- 
verting mode. It is a simple extension to include the 
phase response of the amplifier. 

When the switch in Figure [8] is closed to complete the 
feedback loop, input current is split between the SQUID 
coil and the feedback resistor, the latter of which has 
its impedance reduced by the forward gain of the circuit 
so that the total input impedance is the parallel 
combination 



7SQ\ 
'IN I 



R 



FB 



1 — GpWD 



(18) 



For this system G^^d > across a 1 MHz bandwidth 
and a computer controlled switch allows feedback resis- 
tors RpB = 3.3k57, 5k57, or lOkil to be used. The cir- 
cuit input impedance is 22-69 mD, in the frequency range 
300-900 kHz for the Rpp ^ 10kl7 feedback setting most 
commonly used for SPT and APEX-SZ. 



The SQUID response function (Equation 16) defines 



its open-loop dynamic range as 
The shunt-feedback loop-gain 
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defines the extent to which the current through the 
SQUID coil is canceled by negative feedback, extending 
the SQUID dynamic range to 



rmax 
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loop 



(20) 



and vastly improving its linearity. For the fMUX imple- 
mentation, A^^p ~ 38, 25, 12.5 for Rfb ^ 3.3kn, 5kf7, 
and lOkn respectively. This results in a dynamic range 
that is extended by a factor 23, 15, and 7, respectively. 

When the SQUID dynamic range is exceeded, the cir- 
cuit "fiux-jumps" such that the SQUID is locked in a 
state one fiux quantum or more away from its original 
locking point. This is similar to an operational amplifier 
hitting its rail and 'sticking' there. The FLL has dras- 
tically deteriorated dynamic range and linearity in its 
fiux-jumped state. A fiux-jump can be caused by elec- 
trical pick-up (typically of fast, digital signals), electrical 
glitches, or (very rarely) by a large magnitude electronic 
noise fiuctuation. 

For the 120 SQUID arrays in SPT, the flux-jumping 
rate is typically between to 1 SQUID arrays per 24 hour 
observing cycle. There are several means of recovering 
from a fiux-jump. The simplest method is to 'reset' the 
loop by either opening and closing it again (using the 
switch located in the feedback path, Figure [8| or remov- 
ing and reinstating the SQUID bias. These 'reset' meth- 
ods temporarily set the loop-gain A^^^ to zero. This 're- 
set' only works if the input current through the SQUID 
coil is zero, otherwise the SQUID-loop will not be locked 
at the original locking point. A complication for this 
method is that, with ^j^^p = 0, the input impedance 
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of the FLL is altered, which alters the resonant fre- 
quencies of the LCRhoio resonant circuits. This can spoil 
ETF, causing the bolometer to oscillate and latch in a su- 
perconducting state. For the fMUX system, the SQUID 
input is always non-zero, due to sky-signal and imper- 
fections in setting the nulling amplitude. Consequently, 
instead of 'resetting' the FLL, SQUID-jumps are recov- 
ered from by injecting a corrective current through the 
SQUID flux bias wire, and increasing this current until 
the SQUID jumps back to the original locking point. 




3. Flux Locked Loop Stability 

To maintain stability in the SQUID FLL circuit, phase 
shifts along the feedback loop must provide for negative 
feedback up to the frequency where the loop-gain band- 
width product falls below unity, 38 MHz for Rfb — 
3.3kr2. The room temperature operational amplifier, 
with a bandwidth of 1 MHz, induces a 90° phase shift, 
allowing for an additional margin of 45°. This places 
stringent requirements on the length of the wires con- 
necting the SQUIDs to the room temperature amplifier 
since a 0.66 m round-trip propagation delay amounts to 
45° at 38 MHz. A one-way wiring length of 0.2 m is used 
for APEX-SZ and SPT. In the absence of phase shifts 
induced by strays, this lead length would allow stability. 
In practice, phase shifts arise from other non-idealities, 
such as the parallel resonance formed by the capacitive 
coupling of the MUX inductor to its flux-focusing washer 
or the wiring strays discussed in §111 F| 

To improve the FLL stability, a lead-lag filter is intro- 
duced and connects the output of the SQUID through a 
capacitor and resistor to ground, see Figure [8j At high 
frequency, this RC shunt-filter looks purely resistive and 
acts as a voltage divider with the SQUID output, attenu- 
ating the loop-gain. The filter introduces a phase shift at 
intermediate frequency, when the reactance of the capac- 
itor is similar to the filter resistance. The filter capaci- 
tance is 1 nF and its resistance is set by a 0.2 m Manganin 
wire, about 15 f2. This filter introduces a frequency de- 
pendence on the SQUID transimpedance, attenuating it 
by about 30% at the high end (900 kHz) of the fMUX car- 
rier bandwidth. It also introduces a phase shift inside the 
fMUX carrier bandwidth. These two effects contribute to 
a rise in the noise from the room temperature electronics 



at high bias carrier frequency, discussed in { V C 



4- SQUID Electronics Implementatton 

The SQUID electronics are implemented on two circuit 
boards: (1) a SQUID mounting board that houses the 
SQUIDs at 4 K and provides magnetic shielding, and 
(2) a room temperature SQUID Controller that houses 
the FLL operational amplifiers and SQUID bias current 
generation circuits, and conditions signals as they enter 
or exit the cryostat. 



FIG. 9: (Color online) Photo showing the SQUID mount- 
ing board, housing eight NIST SQUID series arrays and the 
bias resistors for the bolometers. The SQUIDs are mounted 
on top of a Niobium film, visible in the photo. A magnetic 
shield enclosure, fabricated from the high magnetic perme- 
ability material Cryoperm, is shown above the circuit board. 



The SQUID mounting boards, shown in Figure |9j 
are 0.15x0.05 m^ custom printed circuit boards hous- 
ing eight SQUIDs each. SQUIDs are sensitive to changes 
in magnetic field, so that effective magnetic shielding is 
essential. Since the system employs series array SQUIDs, 
it is important to attenuate spatial variations as well as 
time variations in the magnetic field. This is because a 
spatial variation across the array would spoil the coher- 
ence of the individual SQUID elements. 

Each SQUID is mounted with rubber cement on top of 
a lOx 10 mm^ niobium film that is epoxied to the SQUID 
mounting board. Niobium is a Type II superconductor 
and serves to pin magnetic field lines, greatly attenu- 
ating their time variability. The entire circuit board is 
placed inside an enclosure built from a material that has 
a high magnetic permeability at cryogenic temperatures 
(Cryoperm, Vacuumschmelze GmbH, Griiner Weg 37 D- 
63450, Hanau, Germany). This attenuates the amplitude 
of the field lines, reducing both the temporal and spatial 
variability. The enclosure is open on two sides through 
narrow "chimneys" to allow wiring connections to the 
room temperature electronics through an edge connec- 
tor, and to the cold multiplexer components through a 
broadside coupled stripline that is soldered directly to 
the board. 

The 30 mil bolometer bias resistors (Ruas in Figure[2]) 
are located on the SQUID mounting boards. This is a 
compromise, because the 4 K temperature means the bias 
resistors contribute about 2.5 pA/-\/Hz of noise referred 
to the SQUID input. This would be reduced to a neg- 
ligible level had this bias resistor been placed at 0.25 K 
with the sub-kclvin multiplexer components. By keeping 
the resistor at 4 K, we are able to connect to each multi- 
plexer comb with just two wires instead of four, reducing 
the heat load on the sub-kelvin stage. 

The SQUID controllers are 180x170 mm^ custom 
printed circuit boards. One board provides the control 
circuitry for eight SQUIDs. Each SQUID control cir- 
cuit includes a high gain-bandwidth operational ampli- 
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fier with low (1 nV/vHz) input noise and two feedback 
resistors {bkfl and lOkil, Rfb in Figure [s]) that can be 
enabled individually or in parallel via software with an 
analog CMOS switch (adjacent to Rfb in Figure |8| that 
is selected for its low charge injection. A DAC provides 
four software programmable voltages for each SQUID 
that are used to (1) provide a current bias, about 150 fiA, 
through the SQUID, (2) provide a flux bias current, 0- 
25 ^A, to the SQUID input coil that is used to set the 
magnetic field through the SQUID, (3) zero the offset of 
the room temperature op-amp, and (4) provide a current 
to a heater resistor next to the SQUID that can be used 
to heat the SQUID above its superconducting transition. 

Analog bolometer bias carrier signals and nulling sig- 
nals are received from the oscillator/demodulator boards 
through isolation transformers on differential twisted pair 
cables, and attenuated on the SQUID controller board 
before being sent into the cryostat. The FLL output sig- 
nals are transmitted from the SQUID controller boards 
to the oscillator/demodulator boards using differential 
drivers feeding twisted pair cables. 

SQUIDs are very sensitive to digital pickup, so the dig- 
ital circuitry on the SQUID controller boards is isolated 
from the analog circuitry by separating the local grounds 
and coupling the connections from the digital to analog 
circuitry through resistor-capacitor filters. Balanced con- 
nections between the SQUID controller and the demodu- 
lator circuit boards also suppress digital interference from 
the crate. Overall, digital noise pickup is negligible. 



5. SQUID Setup and Tuning 
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Current Through SQUID Coil [\iA] 

FIG. 10: (Color online) The SQUID voltage response to 
a linearly increasing current through the SQUID input coil 
(called a "V-phi" curve) is shown for several choices of the 
SQUID bias current ranging from 80 fiA (bottom curve) to 
160 (top curve). The chosen SQUID operating point is 
marked by a closed circle. 




The sorption fridge that cools the sub-kelvin stage of 
the system has a hold time of about 36 hours. Each time 
this fridge is cycled, the SQUIDs are tuned to configure 
them at their optimum operating point. 

As the mainplate is cooled to 4 K, the SQUIDs cross 
their superconducting transition relatively slowly, with a 
thermal gradient across the series arrays. This gradient 
can transform time varying magnetic fields into spatially 
varying trapped flux across the series array. The first 
step in the tuning process is to raise the temperature 
of each SQUID by providing current to a heater resistor 
located next to it, bringing the device into the normal 
state and then allowing it to cool quickly, reducing the 
probability that time varying fields will be converted into 
spatially trapped flux. 

The next step is to map out the SQUID output volt- 
age response to a current through its input coil (approx- 
imated by Equation 16 1 as a function of the SQUID bias 



current, as shown in Figure 10 The bias current choice 
defines which of these "V-phi" curves the SQUID oper- 
ates on. 

The derivative of the V-phi curve is the SQUID tran- 
simpedance Zf^^^. It defines the SQUID small signal 
response and is the quantity that refers the room temper- 
ature electronics noise back to an equivalent noise cur- 



FIG. 11: Block diagram of the demodulator circuit. 

rent through the bolometer. Ideally, ^i^ans would be 
maximized. The peak-to-peak amplitude relates to the 
squid's large signal response. The distance between 
peaks in each V-phi curve is the current required to pro- 
duce a flux quantum through the SQUID coil, and is 
independent of the bias current. For stable SQUID op- 
eration, a smooth curve with a high level of symmetry 
around the chosen flux bias point is preferred. A typ- 
ical choice for the SQUID operating point is indicated 
by the small red circle on the 150 /iA SQUID bias curve 
(2nd from top). This represents a trade off between the 
criteria listed above. 

Once the SQUID response has been mapped and the 
bias currents chosen, the feedback loop is closed so that 
the FLL is in its low input impedance state, allowing for 
bolometer operation. 



C. Demodulator and Digitization 

Each bolometer channel has its own analog demodu- 
lator that locks into the bias carrier and mixes the sky- 
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signal down to base-band. This demodulator is refer- 
enced with the same analog sine wave that biases its 
bolometer, serving to cancel the clock-jitter sideband 
noise. Sixteen demodulator channels are housed on each 



oscillator/demodulator board, as described in [TV A 



Eight demodulators (half of an oscillator/demodulator 
board) monitor each multiplexer module, consisting of 
seven bolometer channels. Seven of the demodulators are 
phase-locked to the in-phase / component of the seven 
bolometer carriers. The system was originally conceived 
with 8 bolometers per module, but this has been reduced 
to seven so that the eighth "helper" demodulator can 
be used periodically to measure the Q component (90 
degrees out of phase) of each channel. Being able to 
simultaneously observe both / and Q components facili- 
tates fast setup of the nulling signal amplitude and phase 
discussed in §IVE[ 

Each module of eight demodulators is AC coupled to 
the differential SQUID controller output through capac- 
itors. An input transformer suppresses common-mode 
signals. In front of the coupling capacitors, a separate 
analog to digital converter (ADC) samples the SQUID 
static offset at 1 kHz, which is necessary for tuning the 
SQUIDs and monitoring for flux-jumps. 




FIG. 12: (Color online) Photo showing the three readout 
electronics crates (top half of each rack) housing the oscil- 
lator/demodulator circuit boards in the cabin of the South 
Pole Telescope. Control computers, cryogenic housekeeping, 
timing, and network electronics are also visible in the bottom 
half of the racks. 



D. Readout Electronics Crates, Interface, and 
Power Consumption 



The signal from an individual bolometer is mixed down 
to base-band by multiplying it with a unit amplitude 
square wave. This is achieved (Figure 11 ) by passing the 



reference carrier sinusoid through a comparator and us- 
ing its output to control an analog switch that acts as a 
synchronous rectifier on the input waveform. The output 
of this rectifier is sent to a pair of low-pass RC integra- 
tor circuits that remove the high frequency out-of-band 
signals produced by the other bolometer channels. The 
rectifier operates in a doubly-balanced mode, with differ- 
ential input signals, not shown in Figure [TT] as well as 
the differential output. This implementation is built up 
from discrete components and has excellent distortion 
performance because there are no non-linear elements. 
An instrumentation amplifier, with gain that is software 
selectable over two orders of magnitude, takes the differ- 
ence between the two integrators. The preceding stages 
have sufficient gain to override the instrumentation am- 
plifier's low-frequency noise. The amplified base-band 
signal is anti-alias filtered with an 8-pole active low-pass 
filter operating at 400 Hz and digitized with a 12-bit ADC 
at 1 kHz. 

A Field Programmable Gate Array (FPGA) on each 
oscillator/demodulator board assembles the outputs from 
the ADCs and sends the data via an 8-bit parallel low 
voltage differential signal (LVDS) connection to the read- 
out control computer. A software low-pass filter with a 
35 Hz pass-band is applied to the data, which is then 
down-sampled to 100 Hz and written to disk. 



The 16-channel oscillator/demodulator circuit boards 
(Figure 6| are 340 x 360 mm^ in size and 20 circuit boards 
are packaged together in a 9U Versa Module Eurocard 
(VME) crate with custom backplanes that provide power 
to the boards and connect to the control computer. Each 
crate is convectively cooled with fans and has an inte- 
grated power supply. The power consumption is about 
900 W per crate including the power delivered to the 
SQUID Controller electronics described above. Three 
crates (2700 W) house the oscillator/demodulator elec- 
tronics for the 840 bolometer channels of SPT (photo 
Figure [l2|. The 280 channel APEX-SZ system requires 
one crate (900 W). 

The center slot of each crate is used for a clock dis- 
tribution circuit board that sends a central 40 MHz sys- 
tem clock and a once per second GPS-derived timestamp 
across the backplane to each oscillator/demodulator 
board. The system clock and timestamp signals can be 
daisy-chained across several clock distribution boards in 
different crates. 

The system is interfaced to a control computer, run- 
ning the Scientific Linux operating system. Command 
signals are received from the control computer through 
a single RS485 serial signal distributed across the back- 
plane and encoded using the communications protocol 
Modbus (see |http:7/www.modbus.org[ ) . Bolometer data, 
sampled at 1 kHz, flow from the readout crate backplanes 
to the control computer through an 8-bit parallel low 
voltage differential signal (LVDS) digital output bus op- 
erating at 40 MHz. 
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A digital token is passed sequentially from board to 
board through a dedicated hardware port to define which 
circuit board has control of the output bus. This protocol 
allows all boards in the system to share the same digital 
output. 



E. Bolometer Tuning 

The bolometers are tuned by adjusting the bias volt- 
age for each detector and measuring changes in current 
through the bolometer using the oscillator/demodulator 
boards. The tuning algorithms are implemented in soft- 
ware running on the readout control computer. 

The bolometer tuning takes place after the SQUIDs 



have been tuned (see ^IVB5) and the LC Rhoio reso- 



nances have been determined with a network analysis 
(see pnE] ). 

The first step is to cool the bolometer stage to about 
~800 mK, which is above the TES transition temperature 
but below that of the aluminum leads connecting the 
detectors to the cold multiplexer. 

Once this temperature is reached, the detectors are 
provided with a sufficiently large electrical bias that the 
devices will remain above their transition temperature 
when the bolometer stage is cooled to its operating tem- 
perature of 250 mK. Since the SQUID dynamic range is 
not sufhcient to handle seven bias carriers simultaneously, 
the seven detectors in a multiplexer module are biased 
sequentially. After each bias is turned on, a nulling sinu- 
soid is adjusted to cancel the carrier signal at the SQUID 
input, before the next TES bias is activated. 

The optimum nulling sinusoid amplitude and phase 
are determined using measurements of both the / and 
Q components of the demodulated bolometer current as 
provided by each bolometer's demodulator and the one- 
per-multiplexer module "helper" demodulator. After the 
optimum phase is determined, the "helper" demodula- 
tor is no longer needed, and can later be used to tune 
another channel. Once all seven channels in the multi- 
plexer module are biased and nulled, the bolometer stage 
temperature is lowered to its operating temperature of 
^250 mK. The nulling amplitudes are fine-tuned at this 
temperature to correct for slight impedance changes in 
the circuit. 

The seven bolometers in each module are lowered into 
their transition sequentially by decreasing the electrical 
bias power. For a given bolometer, the bias voltage am- 
plitude is decreased while iteratively adjusting the nuller 
amplitude and phase to keep the current through the 
SQUID zero (this corrects for the voltage drop across 
the SQUID input, p5 in Figurejlj as discussed in § IIIFl. 
This traces out the bolometer current vs. voltage curve 
(Figure 13 1. In the upper-right portion of the curve, the 
TES electrical bias power is sufhcient to hold the de- 
vice above its superconducting transition. As the bias 
voltage is lowered, the device falls into its transition and 
is held there by ETF. "Turn-around" is defined as the 
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FIG. 13: A current vs. voltage curve is shown for an 
SPT TES bolometer operating dark at a bath temperature 
of 0.29 K. In the upper-right portion of the curve, the TES 
electrical bias power is sufficient to hold the device above its 
superconducting transition, with constant resistance AV/ AI. 
As the bias voltage is lowered, the device falls into its tran- 
sition and is held there by electro-thermal feedback (ETF). 
"Turn-around" is defined as the minimum in the curve, about 
4.7 fiV, when the dynamic resistance dV/dl transitions from 
positive to negative. 



minimum in the curve (4.7 /zV in Figure 13) when the 
dynamic impedance dV/dl transitions from positive to 
negative (the average TES resistance A]// A/ is always 
positive). At this point the detector ETF loop-gain C is 
unity. Once the TES reaches its target operating resis- 
tance, the bias voltage is held constant and the current is 
once again nulled at the SQUID input. All seven bolome- 
ters in the multiplexer comb are tuned sequentially in 
this manner. The tuning of bolometers in different multi- 
plexer modules across the system is performed in parallel. 
For astronomical observations with SPT (APEX-SZ), the 
TES sensors are typically operated at 60-80% (80-90%) 
of their normal resistance corresponding to typical ETF 
loop-gains £^15-20. 

After the electrical bias power of each detector in the 
system has been decreased so the detectors are operat- 
ing in their superconducting transition, the demodulator 
gain is increased from its high dynamic range setting to 
one of its low noise settings, that has 33-100 times larger 
gain. The nulling amplitudes are fine-tuned and astro- 
nomical observations can commence. 

When detector loading conditions change, for exam- 
ple, when a large change in observing elevation is made, 
adjustments in the nulling may be necessary. Significant 
changes in detector loading may necessitate re-tuning of 
the detectors. 
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SYSTEM PERFORMANCE 



B. Electrical Cross- Talk 



In tiiis section, measurements of the system 
performance — including noise, sensitivity, and cross- 
talk — are presented. 



A. System Setup 



For SPT and APEX- SZ, th e bolometer tuning (de- 
scribed in more detail in ' 



IV E ) is integrated with the cy- 



cling of the sub-kelvin sorption fridge to minimize the to- 
tal system setup time. The sorption fridge is cycled every 
24-36 hours, elevating the temperature of the bolometers 
and SQUIDs to above their superconducting transition. 
This necessitates a re-tuning of the bias parameters for 
both. The process used for SPT is outlined below. 

The SQUIDs are returned to their operating tempera- 
ture mid-way through the cryogenic cycle and are tuned 



(§ IVB5) in parallel with the second half of the cycle. 



At this point, the bolometers are above their supercon- 
ducting transition. The SQUIDs are first heated for 30 
seconds to release trapped flux and are subsequently al- 
lowed to cool quickly through their transition. Following 
this cooling, it takes about 25 minutes for the SQUIDs 
to return to their base operating temperature. Next, 
voltage-current curves (Figure 10 ) are measured for each 



SQUID to determine their optimal flux and current bias. 
For SPT, tuning aU 120 SQUIDs takes ~25 minutes, and 
is performed in parallel with the sorption fridge's cycle, 
therefore coming at no cost to the observing efficiency. 

Towards the end of the cryogenic cycle, the bolometer 
stage is left at ^800 mK so that the detectors are above 
their superconducting transition. The control computer 
commands the readout system to bias the detectors with 
enough electrical power to hold the TES normal and tune 
the nulling current to cancel the bias at the SQUID in- 
put (j jlVE ); this takes about 10 minutes. After this, it 
takes about 2 hours for the sorption fridge to cool the 
bolometer stage to its ^250 mK base temperature. Once 
the cooling is complete, it takes about 15 minutes to 
adjust each detector voltage bias to the desired point 
in its superconducting transition, and then to re-adjust 
the nulling current (j |IVE |. An I-V curve is recorded for 
each detector (e.g.. Figure 13 1. The total time elapsed 



for readout /bolometer system tuning and setup is about 
25 minutes, not counting the time (several hours) that 
would otherwise be spent cycling the cryogenics. 

SPT has 724 operational detectors. Others have wiring 
or other defects that prevent their operation. Typically 
650±15 of these detectors show good performance each 
tuning and are used for astronomical observations. De- 
tailed statistics and characterization of end-to-end re- 
ceiver performance are presented in a separate publica- 
tion i52l. 



Channel-to-channel cross-talk can occur through op- 
tical signals or electrical signals. Optical cross-talk be- 
tween adjacent bolometers in SPT and APEX-SZ is typi- 
cally about 1% and is determined by the optical design of 
the experiment and the arrangement of detectors within 
a silicon wafer. The goal is for electrical cross-talk aris- 



ing in the readout system and discussed in § III C to be 
small in comparison to this optical cross-talk. 

Electrical cross-talk in the multiplexer readout has 
been measured both in the laboratory and during sky- 
observations. For the laboratory measurement (25] . 
one detector, mounted in its own light-tight enclosure, 
was excited with light from an LED. Five other sensors 
were read out on the same fMUX module, but mounted 
in a separate light-tight enclosure. The absence of any 
measurable signal in the five dark channels established 
an upper limit of 0.4% electrical cross-talk. 

Cross-talk has been characterized with on-sky obser- 
vations of RCW38, MAT5A, Venus and Jupiter with the 
South Pole Telescope. Though we cannot distinguish be- 
tween positive cross-talk that arises from optical versus 
electrical sources, negative cross-talk should be domi- 
nated by electrical signals. The signal-to-noise levels for 
RCW38 and MAT5A are substantially lower than for the 
planets, and no negative cross-talk is observed. Using 
the observations of Venus and Jupiter, which produce 
signals that are large enough for small cross-talk signals 
to be observed in neighbors, a median cross-talk between 
neighbors of 0.3% is observed. This is consistent with 
the expectation of 0.25% and 0.3% from mechanisms (2) 
and (3) described in § III C 



C. Noise Performance 

In this section, the measured noise levels for a dark 
detector in the laboratory and an SPT detector observing 
the sky are compared with the theoretical expectations. 

The noise contributions are divided into three cate- 
gories: photon noise, bolometer noise, and readout sys- 
tem noise. Within each category, a contribution may 
arise as a current or power source. Sky signals arise as 
power sources (photons converted to heat), so that noise 
equivalent power (NEP) is the relevant metric. Other 
sources, such as the SQUID noise, are constant in terms 
of an input referred current, and are labeled here as cur- 
rent noise sources. The NEP of current noise contrib- 
utors depends on the device parameters and operating 
conditions and is converted from current to NEP by mul- 
tiplying the SQUID-input referred current noise by the 
inverse of the detector responsivity, ~ ~^hSs^ / (see 



Equation A6 in Appendix Ah which is valid in the limit 



of high loop-gain. 

Achieving optimal noise performance with TES 
bolometers for astronomical observations requires a care- 
ful optimization of detector and readout system param- 
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eters. Most important amongst these is the power re- 
quired to saturate the detectors, -Psat- If this power is 
many times higher than the absorbed radiation power (or 
the detector is dark, so that there is no radiation power) , 
thermal phonon noise will be high and relatively large 
voltage biases must be used to provide the Joule heat- 
ing necessary to keep the detector in its transition. The 
high voltage bias results in an amplification of current 
noise sources when they are referred to NEP, so that the 
readout component of the noise is greatly enhanced. On 
the other hand, if Pgat is too small, the detector will not 
have enough dynamic range to operate through a broad 
range of observing conditions. Typically Psat is set to be 
roughly twice the expected absorbed incident radiation 

(« 2 X P,ad) ■ 

Before describing the noise contributions below, we 
note that a mixer responds differently to noise that is 
encoded as a modulation on the carrier and noise that is 
superimposed on the carrier at the same frequency (see 
Appendix |A]) . This implemen t at ion of the fMUX system 
uses a square- wave mixer (§ IV C I, which introduces a 



further distinction between (a) broadband noise sources 
that are superimposed on both the carrier frequency and 
its harmonics, and (b) narrow-band noise sources that 
are superimposed on the carrier but do not extend out 
to the carrier harmonics. SQUID noise is an example 
of (a) because its bandwidth is not limited. Bolometer 
Johnson noise is an example of (b) because its bandwidth 
is limited by the LCRhoio filter. For a square mixer, 
these noise sources are enhanced in the demodulation 
process by factors of | and respectively, in compar- 
ison to noise that modulates the carrier. Fortunately, 
in an alternating-voltage biased system, power terms re- 
ceive an enhancement by a factor y/2 (with respect to 
a constant-voltage biased system) during the modula- 
tion process, so that the signal-to-noise ratio is nearly 
unaffected by these factors in comparison to a constant- 
voltage biased system. Appendix [A| presents a detailed 
discussion and derivation of the relevant factors. 



Dark Notse 



In Figure 14 (top, solid line) the noise equivalent power 
(NEP) density is shown for a representative SPT detec- 
tor operated in the laboratory with no incident radiation 
power ("dark"). The detector has a saturation power of 
GAT ~ 25 pW. Its normal resistance -RNormai = 0.98 
and it has been lowered into the transition to 0.73 i?Normai 
with a voltage bias of 4.2 /iVrms (the current- voltage 
curve for this detector was shown in Figure 13). The 



roll-off in the measured noise at 35 Hz is due to the read 
out system's software low-pass filter described in 



IV C 



Noise expectations for this detector are summarized in 
Table |llj Columns are included for dark operation (zero 
radiation power, shown in the top panel of Figure [M] 
as a dashed line) and for a radiation loading of 11 pW 
(shown in the bottom panel of Figure 14 and discussed 
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FIG. 14: The measured noise equivalent power (NEP) 
spectrum density for an SPT TES bolometer operated at 
0.73i?Normai in the laboratory with no incident radiation 
( "dark" ) is shown (top, solid line) with the theoretical expec- 
tation from Table superimposed (top, dashed line). This 
configuration represents the worst possible configuration for 
the readout noise contribution, since the incident power is 
zero. This means the voltage bias, which refers current noise 
sources to noise equivalent power, is maximal. In the lower 
panel the on-sky noise for this detector is estimated assuming 
a total radiation loading of 11 pW and the detector param- 
eters of Table |ll] With radiation power, the detector voltage 
bias is lower and the readout noise is reduced when referred 
to NEP. 



V C 2[ ), which is typical of the 150 GHz detectors for 



SPT. Table [H] also includes the equations and parameters 
used to estimate each noise component. A similar table 
appears in Ref. [2] for APEX-SZ detectors. The noise 
contributions for "dark" operation are described below. 

TES detectors exhibit thermal carrier (or phonon) 
noise, which arises from random variations in the flow 
of phonons from the TES to the thermal bath. A factor 
7ne has been included in the Table [H] equation describ- 
ing this noise to account for a temperature gradient along 
the thermal link |M]. We have taken 7ne — 0.5 56J, con- 
sistent with a normal metal link at our operating tem- 
peratures. For the dark detector power spectrum den- 
sity shown in Figure [M} thermal carrier noise dominates 
above 1 Hz and is eventually cut-off by the TES time 
constant t^es ~ 1 ms, which is not visible in the plot as 
it is beyond the low-pass filter cut-off. 

Both the detector and bias resistor produce Johnson 
noise in the narrow LCRhoio pass-band that surrounds 
the carrier. The TES detector Johnson noise current 
■y/ 4/cBTboio/-Rboio is superimposed on the amplitude mod- 
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Detector Parameter: 

TES Transition Temperature rboio=550mK 
Thermal Conductivity 

Center Frequency & Bandwidth =153 GHz 

Detector Setup: 

Bath Temperature 
Incident Radiation Power 

TES Normal and operating Resistance -RNormai = 0.98 f2 
Bias Voltage 
Loop-Gain 

Noise Source Equation 



G = If ~ 96pW/K, gT = ^~ 165pW/K 
Au =38 GHz 

Dark 
Tbath=290mK 

Prad =OpW 

^boio ~ 0.73 

Vbias=4. 2 /iVflMS 

£~ 25 



On-Sky 
rbath=260mK 
P,ad = llpW 

-Rboio = 0.73 
Vbias=3.2 /iVrms 



£~ 15 



NEP (aW/VH^) 



Photons: 

Shot noise 
Correlation noise 



\/2/li/Prad 



Dark 






On-Sky 

47 
^/?• 56 



Bolometers: 

Thermal Carrier noise (7ne = 0.5) 
Johnson noise 



\/7NE4fcBr^olo5T 

olo-fclcc 



37 

28//; 



37 

21/jC 



Readout: 

Warm electronics 
SQUID 

Bias resistor Johnson noise 



/Hz •/2 
pA Vbia. 
/Hz V2 



3.7- 
3.5 



17 
16 
15 



13 
12 
11 



Total 



47 



64,71,85 
= 0,C = 0.3,C = 1) 



TABLE II: Theoretical noise equivalent power expectation for a typical SPT detector. The "dark" column corresponds to 
the operating conditions for the noise spectra recorded in the laboratory and shown in Figure [M] wherein there is no incident 
radiation. The "on-sky" column corresponds to typical incident radiation power for a 150 GHz SPT detector. Totals are 
included for several cases of photon correlations (^=0, 0.3, and 1). A similar table appears in Ref. [2] for APEX-SZ detectors. 
Device and experiment related parameters are defined at the top of the table, h is Planck's constant, fee is the Boltzmann 
constant, ^ is the photon correlation parameter, and 7ne is the non-equilibrium parameter [64j . 



ulated carrier. For a sinusoidal-biased TES, this Johnson 
noise beats with the carrier and is mixed down to base- 
band where its thermal power lies within the detector 
time constant and is suppressed by electro-thermal feed- 
back. This suppression is similar to that which takes 
place for a static-biased TES [61], and is derived in Ap- 
pendix B of Ref. [55]. For the operating parameters of 
SPT and APEX-SZ detectors, the TES Johnson noise 
suppression needs to be calculated numerically, but can 
be approximated by ~ 1/C for the component that is 
in-phase with the bias carrier |65| (the out-of-phase com- 
ponent receives no suppression, but is discarded by the 
demodulator). The bias resistor suppression is less than 
20% [5S] and has been neglected in the noise expecta- 
tion tabulated in Table HT] For a sinusoidal-biased de- 
tector, Johnson noise is enhanced by a factor be- 
cause uncorrelated noise from both sidebands appear 



post-demodulation, as discussed in Appendix [Xj This 
noise is referred to a noise equivalent power by the in- 
verted responsivity ~ ~^uas^ The result is that 
the total Johnson noise is V 4/cB7boio-Peioc/'C, exactly as 
it would be for a detector biased with a static voltage. 

The readout electronics contribute noise from the bias 
resistor, SQUID, and room temperature electronics as 
described below. As described in § |IV[ each multiplexer 
module bias resistor is located at 4 K adjacent to the 
SQUIDs (see Figure [T]) rather than on the sub-kelvin 
stage adjacent to the detectors. This higher temperature 
results in higher Johnson noise from the bias resistor. 
The advantage is that the number of wires connecting 
to the sub-kelvin stage is just two per multiplexer comb, 
rather than the four wires that would be required if the 
bias resistors were located on the sub-kelvin stage. This 
noise contribution could be made negligible by using a ca- 
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pacitive voltage divider [SU] , an inductive bias, or moving 
the bias resistor to lower temperature. 

The SQUIDs contribute a white noise level of 
3.5pA/VH^. Since the detectors are biased at high fre- 
quency, the bolometer signals are modulated at frequen- 
cies well above the 1 /f noise of the SQUIDs and pre-mixer 
amplifiers in the room temperature system. The room 
temperature electronics noise is 3.7pA/-\/Hz, after being 
referred back to the bolometer current by the SQUID 
transimpedance, zf^^^. A higher ^t^ans would reduce 



this noise that is dominated by the 1 nV/ vHz input noise 
of the first stage amplifier that follows the SQUID inside 
the flux locked loop. The room temperature electron- 
ics noise also receives contributions from Johnson noise 
in the resistors that set the gain of the amplifiers and 
convert voltages to currents. 

When we account for these noise contributions, the 
total noise expectation is in good agreement with mea- 
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surements of dark detectors as shown in Figure 14 



FIG. 15: Power spectrum distribution of a bolometer ob- 
serving in the 150 GHz band on the South Pole Telescope is 
shown. The low frequency noise is primarily due to atmo- 
spheric fluctuations. This detector's performance is typical 
for detectors in the 150 GHz band and agrees well with the 
expectation, shown as a shaded band. Its measured noise 
equivalent temperature (NETcmb) is 409 /xK • projected 
on the sky. 



2. On-Sky Noise 

The noise expectation for the detector shown in Fig- 
ure [M] extrapolated to typical SPT 150 GHz observing 
conditions with llpW of incident radiation power, is 
shown in the "on-sky" column of Table [n] This radia- 
tion power results in a reduction of the Joule heating 
necessary to keep the detector at 0.73i?Normah so that 
the detector is tuned with a lower bias voltage. This 
lower bias voltage results in a reduced NEP for all current 
noise sources, such as the three components of readout 
noise listed in Table |Tl] Since the low frequency noise 
from the readout system is proportional to the bias volt- 
age amplitude, its amplitude is also reduced. If readout 
were the only contribution to low-frequency noise, the 
knee would be at 0.4 Hz for this detector parameters. 

In addition to the thermal carrier noise, ETF sup- 
pressed Johnson noise, and readout noise, fluctuations 
in the rate of incident photons produces an additional 
noise term corresponding to photon shot noise and pho- 
ton correlations [66, 67 , also included in Tabic Iljfor the 
"on-sky" detector. The degree of photon correlation de- 
pends on poorly characterized details of the filters and 
coupling, so that the coefficient ^ [55] is often introduced 
to parametrize a correction factor for the simple equation 
presented in Table [H] ^ lies between 0-1, with 0.3 being 
a reasonable choice. Table [Tl| includes noise expectations 
spanning the full range. 

The power spectrum distribution for a representative 
150 GHz bolometer operating "on-sky" during the 2009 
observing season on the South Pole Telescope is shown in 
Figure [15] This data corresponds to the raw timestream; 
no filtering has been applied to remove atmospheric sig- 
nals or other effects. The SPT bolometer stage was op- 
erated at a bath temperature of Tbath =260 mK and the 
detector has sufficiently similar device parameters to the 



dark detector, shown in Figure 14 that its total noise ex- 



pectation is the same to within 2% of the values presented 
in Table [T4[ Superimposed on Figure [T5| is a shaded band 
corresponding to the noise expectation for this detector, 
excluding atmospheric contributions. The width of the 
shaded band indicates the range of allowed photon cor- 
relations < ^ < 1. 

At frequencies below ^3 Hz, the spectrum is domi- 
nated by atmospheric fluctuations that arc not included 
in the theoretical expectation presented in Table [Tl] For 
science analyses, filtering and subtraction techniques (see 
e.g., [33]) applied in the data processing pipeline are 
used to reduce the effect of atmospheric signals to be- 
low 1 Hz. Above this frequency, the spectrum agrees 
well with the expectation. After calibrating this detector 
response with astrophysical sources, its measured sensi- 
tivity is 409 /iKcMB ■ -v/s projected on the sky and referred 
to a source at the CMB temperature. This is typical of 
SPT 150 GHz detector performance. Noise and sensitiv- 
ity statistics for the APEX-SZ and SPT receiver systems 
are presented in Refs. (52J and [2^ respectively. 



VI. CONCLUSIONS 

We have presented the design and performance of the 
SQUID-based frequency-domain multiplexed (fMUX) 
readout system that has been deployed for TES detectors 
on the APEX-SZ instrument and the South Pole Tele- 
scope SZ-camera. The readout reduces the heat-load on 
the sub-kelvin cryogenic stage by sharing wiring among 
many detectors using a simple architecture that mini- 
mizes the sub-kelvin system complexity. It has near-zero 
power dissipation on the sub-kelvin stage, is extensible 
to much higher multiplexing factors, allows bias levels 
to be set individually for each detector, is modular with 
no shared components between readout modules, and ex- 
hibits strong reduction of microphonic response and low- 
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frequency EMI because the sky signals are modulated up 
to high frequency. 

The system has been characterized in terms of noise 
performance, sensitivity, detector yield, and cross-talk in 
both the laboratory and with on-sky observations. The 
measured noise equivalent power of the system is con- 
sistent with theoretical expectations, and contributions 
from the readout system are small in comparison to de- 
tector and photon noise above 1 Hz. At low frequencies, 
the system noise is dominated by residual atmospheric 
fluctuations. 

These implementations have shown that frequency- 
domain multiplexing is an effective technology for scal- 
ing to very large arrays of TES detectors. An upgraded 
digital version of the analog fMUX back-end electron- 
ics has recently been developed [5Tj and has been de- 
ployed for the EBEX [55] balloon-borne CMB polarime- 
ter and the ground-based POLARBEAR ^ CMB po- 
larimeter. A polarimeter camera |71| for the South 
Pole Telescope is being commissioned during the austral 
2011/2012 summer and uses the digital fMUX. The 
digital fMUX system achieves an order of magnitude 
lower power consumption, higher multiplexing factors, 
faster detector and SQUID setup, better low frequency 
noise performance, and lower cost. Development is un- 
derway to develop fMUX electronics for satellite appli- 
cations [75]. Present prototypes use Digital Active 
Nulling [73^ in place of the flux-locked loop. This extends 
the SQUID bandwidth to allow for 64 bolometers per 
SQUID module using a single pair of sub-Kelvin wires. 



Acknowledgments 



We thank Andy Smith for useful discussions and de- 
sign/fabrication of niobium inductors at Northrup Grum- 
man, Kent Irwin and Gene Hilton for useful discussions 
and design/fabrication of the series array SQUIDs at 
NIST, and Hannes Hubmayr, Bryan Steinbach, Kam 
Arnold, and Graeme Smecher for useful discussions of 
bolometer responsivity and demodulator response. 

The National Science Foundation funds APEX-SZ 
through grants AST-0138348 & AST-0709497 and the 
South Pole Telescope through grants ANT-0638937 and 
ANT-0130612. Work at LBNL is supported by the 
U.S. Department of Energy under Gontract No. DE- 
AC02-05CH11231. The McGill team acknowledges fund- 
ing from the Natural Sciences and Engineering Research 
Council of Canada, Canadian Institute for Advanced Re- 
search, and the Canadian Foundation for Innovation. 
MD acknowledges support from the Canada Research 
Chairs program and a Sloan Fellowship. 



Appendix A: Modulation and demodulation transfer 
functions for alternating-voltage biased bolometers 

In this section, we consider the transfer functions for 
the modulation and demodulation of signals for a bolome- 
ter voltage biased with a constant amplitude sinusoidal 
carrier. We also consider the transfer function for noise 
that is superimposed on the amplitude-modulated car- 
rier. 

First consider a single Fourier component of a small 
power signal PsCosiOst absorbed by a bolometer that is 
biased deep in its transition with a constant amplitude 
voltage bias Vbias = Vc sin ojct. 

Under the assumption of strong electro-thermal feed- 
back (ETF) , the total power on the detector will be con- 
stant and determined by the average thermal conductiv- 
ity G and temperature difference from the bolometer to 
its bath AT = Tboio - Tbath, 

frOTAL = GAT = Vbias/bolo + Prad + Ps COSCJ.i, (Al) 

where Trad is the approximately constant radiation load- 
ing power from the sky. The amplitude of the bolometer 
voltage bias is constant, so that the power signal pro- 
duces an amplitude modulation Ig on the bolometer car- 
rier current T: 



(A2) 



4oio = {Ic + Is COS uJst) sin ujct- 



Keeping only the time varying terms from Equation |Al[ 
the power balance is 

— Ps COSUJst = Vcis COS Wst SIT? [iUct). (A-3) 

Since the carrier frequency Wc is much faster than the 
bolometer time constant, the sin^(wcO factor integrates 
to 1/2, leading to 



(A4) 



Thus, the strong ETF responsivity for a bolometer biased 
with a sinusoidal voltage is 



V2 

^RMS 



(A5) 



for amplitude modulated signals. For simplicity, we 
have left out the factor that accounts for the bolome- 
ter time constant. Equation |A5| should be compared to 
^ ~ —y- for a constant- voltage biased bolometer [71] . 
For finite loop-gain C, the responsivity in Equation |A5| 
becomes 



dPs 



V2 



yRMS 



(A6) 



Next we consider a bias carrier waveform Ic sin ujct 
and present the response functions for sine-wave demod- 
ulation (the traditional mixer reference waveform) and 
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FIG. 16: A generic demodulator, consisting of a mixer and 
low-pass filter is shown (left). Waveforms for the special case 
of a sinusoidal input demodulated with a square wave refer- 
ence is shown at right. 



square-wave demodulation (the mixer reference used for 
this implementation of the fMUX system) as a function of 
the phase difference between the carrier and the demod- 
ulator, (/)rct. The results for in-phase / demodulation and 
90-degree out-of-phase Q demodulation are summarized 
in Table HE 



A generic demodulator is shown in Figure 16 (left), 
where an input signal x{t) is multiplied by a reference 
signal c{t) then integrated with a low-pass filter (LPF) to 
produce the output z{t). The cutoff frequency of the LPF 
is much lower than the reference frequency and much 
larger than the modulation-encoded sky-signals. 

For the special case of an on-resonance, in-phase sinu- 
soidal input x{t) = Ic sm.{ujct) and a unit amplitude sine- 
wave reference c(t) = sin(a;ct), the mixer output is the 
familiar y{t) = x(t) x c{t) = hsin^{ujct) = 
The LPF averages the second term to zero, so that the 
demodulator baseband output is z{t) — Ic/^. 

The x{t), c{t), y{t), and z{t) waveforms are shown 
for this special case of a square-wave mixer in Figure 16 
(right). The reference is an on-resonance, in-phase, unit 
amplitude square wave c(t). The output is a baseband 
signal with amplitude equal to the average of a rectified 
sine wave, ^Ic- If the phase of the reference differs from 
the sinusoid phase by ^rof, the amplitude of the base- 
band output signal z(t) = ^/^ cos (/iicf , meaning that the 
response function is - for in-phase / demodulation and 
zero for Q. 

Figure [T7| shows other demodulator input signals that 
are relevant for understanding the signal and noise trans- 
fer functions of the fMUX system. The first case (a) 
corresponds to a bias carrier Ic sinwc^, 3> wlpf mod- 
ulated by a small sky signal Ig cosw^i, ujg <^ i^lpFj yield- 
ing a demodulator input signal 



x{t) — (Ic + Is coswst) sin Wei 



(A7) 



= Ic s'mujct + [sin(a;c -|- UJg) + sin(a;c — uJs)] ■ 

This corresponds to the amplitude modulation produced 
by a power signal discussed in the first part of this Ap- 



FIG. 17: Frequency-domain representations of (a) a modula- 
tion of carrier amplitude Ic with a signal of amplitude Is , (b) 
an uncorrelated sinusoid with amplitude /test adjacent to car- 
rier amplitude Ic, (c) narrow-band white noise with amplitude 
*noise superimposed at frequencies surrounding carrier ampli- 
tude Ic, and (d) broadband white noise with amplitude jjjoise 
superimposed on carrier amplitude Ic- A fifth case, where 
narrow-band white noise modulates the carrier (as would be 
the case for bolometer phonon or photon noise) is not shown. 



pendix; and arises from sky signals, bolometer phonon 
noise, and bolometer photon noise. The sky signal is 
split into correlated sidebands adjacent to the carrier and 
appears in the demodulated timestream with amplitude 

z{t) = ^ Is cos cos UJgt. 



Case (b) in Figure 17 corresponds to an off-resonance 
sinusoid It sin ujtt injected adjacent to the carrier. For 
simplicity, this test sinusoid is assumed to be injected af- 
ter the bolometer so that the sinusoid does not deposit 
power in the bolometer that would be subject to ETF. 
The test sinusoid beats with the carrier so that the de- 
modulated signal is z{t) — ^/t cos 0rcf sin |wt — ujc\t for 
a sine- wave mixer and z{t) = ^/t cos 0rof sin jwt — ujc\t 
for a square-wave mixer. Notice that the test signal 
has the same demodulation attenuation factor as an on- 
resonance amplitude modulated signal. Test signals of 
this type are used extensively in the fMUX system for 
measuring electronics transfer functions and characteris- 
ing detectors. 

The situation for white noise added to the carrier wave- 
form (i.e., does not amplitude-modulate the carrier) re- 
quires special attention, because it differs from noise 
that amplitude modulates the carrier and because the 
square- wave mixer responds differently to broadband and 
narrow-band noise sources. 



Case (c) in Figure 17 corresponds to bandwidth lim- 
ited white noise adjacent to the carrier with amplitude 
x{t) — i„oise[A/\/Hz]. An example would be Johnson 
noise from a resistor in an LCR notch filter (or from 
a normal bolometer in our fMUX circuit, though ETF 
complicates the situation for a bolometer in transition, 
suppressing this noise) . A single sideband can be decom- 
posed into sinusoids, and hence has the same response 
function as case (b). There are, however, contributions 
from the sidebands on either side of the carrier, each 
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Input Signal x{t) 



c{i) =Square-wave reference 
/ Q 



c{t) — sm{ujct) reference 
/ Q 



Sinusoidal carrier 

Amplitude modulated carrier (a) 

Sinusoidal test signal (b) 

Superimposed narrowband white noise (c) 

Superimposed broadband white noise (d) 



Ic sin ij ct 
{Ic + I3 cos Ldst) sinuict 
It sinojtt, u!t / 
n{t) 
n{t) 



|/c 

TV 

\'^It sin{ut - LOc)t 
n{t) 



^C^^s COS CJst 

2 



t|/t sin(a;t — ujc)t 



TABLE III: Demodulator response factors are tabulated for various signal and noise types including those shown in panels 
(a) through (d) of Figure 17 uj^ is the demodulator reference frequency. Notice that white noise that is superimposed on the 
carrier (not amplitude modulated) has a factor y/2 enhancement over amplitude modulated signals. For the case of a square 
mixer and broadband white noise, the factor is ^ because noise located within the LPF bandwidth of the reference signal's odd 
harmonics also contribute. 

f An arbitrary phase which depends on the start time is omitted from these equations. 



contributing the same power. The two sidebands are un- 
correlated with one another, resulting in a V2 enhance- 
ment, for a total contribution z{t) = ^inoiso for a sine- 
wave mixer and z{t) = ^^inoisc for a square-wave mixer. 
There is no dependence on ^rcf since this input signal is 
not coherent. The end result is that the sky-signal to 
narrow-band noise ratio is the same for a TES read out 
through a square- wave mixer and a sine- wave mixer, and 
equivalent to what would be obtained for a TES with a 
static bias. 



Case (d) in Figure 17 corresponds to broadband white 



noise with amplitude Znoiso[A/vHz]. An example would 
be SQUID noise or readout electronics noise. For a sine- 
wave mixer, the response is identical to that of case (c). 
For a square-wave mixer, both the noise power adjacent 
to the carrier and the noise power adjacent to the odd- 
harmonics of the carrier contribute. This is immediately 
evident by examining the frequency-space representation 
of the square-wave, with includes delta functions at all of 
the odd- harmonics that fall in amplitude as 1/n. Since 



the noise waveform is characterized by random phase and 
the square-wave mixer simply rotates the phase of half 
the time samples by 180°, it is easy to see that the de- 
modulated amplitude will be unchanged, z(t) = inoisc- 
We note that the demodulated amplitude of broadband 
noise can be reduced to that of case (c) if the carriers 
span less than 2 octaves by placing a low-pass filter in 
front of the mixer, removing input from the 3'"'' harmon- 
ics. The end result is that the sky-signal to narrow-band 
noise ratio is the same for a a TES readout out through a 
square- wave mixer and a sine- wave mixer, and equivalent 
to what would be obtained for a TES with a static bias. 
The end result is that the sky-signal to broadband noise 
ratio is (2/7r)/(l/^) ~ 0.90 worse for a TES read out 
through a square-wave mixer as compared to a sine- wave 
mixer. The latter is equivalent to what would be ob- 
tained for a TES with a static bias. Typically for fMUX 
readout, the contribution from broadband noise sources 
is small, and this difference is negligible. 
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